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Abstract 



< 

' In this work we are concerned with the problem of achieving max-min fairness in Gaussian 

parallel channels with respect to a general performance function, including channel capacity or 
decoding reliability as special cases. As our central results, we characterize the laws which determine 
the value of the achievable max-min fair performance as a function of channel sharing policy and 
^ • power allocation (to channels and users). In particular, we show that the max-min fair performance 

Q , behaves as a specialized version of the Lovasz function, or Delsarte bound, of a certain graph 

induced by channel sharing combinatorics. We also prove that, in addition to such graph, merely 
a certain 2-norm distance dependent on the allowable power allocations and used performance 
^ ■ functions, is sufficient for the characterization of max-min fair performance up to some candidate 

! interval. Our results show also a specific role played by odd cycles in the graph induced by the 

channel sharing policy and we present an interesting relation between max-min fairness in parallel 
channels and optimal throughput in an associated interference channel. 



Index Terms 

Lovasz function, Delsarte bound, parallel channels, max-min fairness, graphs 

I. Introduction 



H I Fairness represents an important goal in the design of power, bandwidth and time allocation 

policies for multi-user channels. It is usually desired to achieve fairness with respect to com- 
munications and information theory metrics, such as spectral efficiency, decoder reliability, 
etc. [1], [2]. The mostly used notion of fairness is the max-min fairness, which is an instance 
of equity in terms of economy markets and consists in the maximal possible improvement 
of the worst performance metric [3] [4], [5]. 

The single-user communication over parallel channels is a well-studied topic both from the 
viewpoint of information theoretic optimality as well as suboptimal practical power allocation 
approaches, see e.g. [6], [7], [8], [9], [10] and references therein. However, the max-min fair 
allocation of power, bandwidth and time to multiple users sharing the parallel channels access 
still poses practical problems and needs a deeper understanding [1], [11]. The issue of max- 
min fairness in (multi-user) parallel channels has been addressed in [11], [12], [13], [14] and 
references therein. Concurrently, a cellular downlink or uplink using Orthogonal Frequency 
Division Multiplex (OFDM) appears to be the most relevant example of parallel channels 
shared among multiple users. In [11], the max-min fair carrier and antenna assignment is 
studied for a multiple antenna OFDM downlink. More generally, in [12], [13] the performance 



of max-min fair power allocation and max-min fair assignment of parallel channels is 
characterized within the framework of so-called blocking and antiblocking polyhedra [15], 
[16]. The characterization of user performance achieved under max-min fairness is provided 
in [12], [13] in the form of bounds and duality-like optimization problems. 

Several important aspect of the fairness problem in parallel channels, such as e.g. the 
optimum interrelations between the combinatorics of channel sharing and the real-valued 
power allocation, are still open in the general case. Also the essential straight questions such 
as "what is the user performance under max-min fairness in parallel channels equal to?" or 
"what is the power/time/bandwidth function describing it?" remain unanswered so far. In 
this work we make a step towards satisfying answers to the above questions in Gaussian 
parallel channels, when the interest is in max-min fairness with respect to user performance 
measured by a general performance/QoS function; this includes the most celebrated cases 
of channel capacity, spectral efficiency, decoder reliability (unity minus decoder error rate), 
etc. We state insightful optimistic and pessimistic bounds on the user performance (Sections 
ITVl rvj ). The essence of our results is that, under constraints on transmit power, the max- 
min fair performance behaves as a specialized version of the Lovasz function [17], [18] of a 
specific graph, which is induced by the channel sharing topology, or combinatorics. We prove 
further that, in addition to such graph description, a certain 2-norm distance dependent on 
the allowable power allocations and users' performance functions is sufficient for enclosing 
the max-min fair performance by some lower and upper bounds. These bounds prove to be 
especially insightful as they offer a separation of influences of the channel sharing policy 
and the real-valued problem of power allocation. We aid the interpretations of the introduced 
channel sharing topologies and the proposed bounds by several parallel channel examples 
and visualizations. 

Our results exhibit a specific role of odd cycles in the graph induced by the channel sharing 
policy. We present also an interesting relation between max-min fairness in parallel channels 
and optimal throughput in an associated interference channel. Furthermore, the presented 
bounds show a link between the user performance under max-min fairness and (zero-error) 
graph capacity [19], [17]. The proofs of our results are constructive and allow for the design of 
several novel power and time allocation algorithms for parallel channels with predetermined 
channel sharing topology; this can be motivated by regulations on processing of traffic 
classes or standardization/hardware constraints (Section rvil). The proposed algorithms offer 
a better performance-complexity trade off than conventional solution methods and achieve 
user performance within some specified distance from the max-min fairness. 

IL Model and Preliminaries 

We consider the set of Gaussian (in the sense of Additive- White-Gaussian-Noise) parallel 
channels, treated as one multi-user channefl. The transmitter-receiver pairs communicating 
with each other over this channel are referred to abstractly as users and are grouped in the set 
/C = {1, . . . , K}. The parallel channels are assumed to be deterministic and frequency-fiat. 

'Notation: The nonnegative and positive orthants in R^^^, wliere we set R^^^^ = R-*^ and C*^'*^ — , are denoted 
as R* *^^ and R;^^^, respectively. By we denote tlie set of symmetric matrices in R^^^ and tiie cones of doubly 
nonnegative matrices and completely positive matrices in R*'''^ are denoted as and P^, respectively (see Appendix 
|A]for the definitions). By St(X) we denote a ball with radius e centered at X € C^^'^ . By >; we denote the usual partial 
order on the set of symmetric matrices and X oY, denotes the Kronecker product of X,Y G C^^^. For any vector 
X = (a;i, . . . , xk) € we define (a;);, — Xk. Given a matrix X £ C^^^ with elements Xki, l<k<K,l<l<N, 
we write simply X = (xki) and define similarly (-X")^.; = Xki- By X' we denote the conjugate transpose of X G C^^^ . 
Given X G C"^^'^, diag{X) G C'^''"^ is such that {diag{X))^^ = [X)^^ and (difflff(X))j^, = Q, k I, 1 < k,l < K . 
Further, given x = (si, . . . ,xk) G R+, a vector X2 is defined as (x^)j, — ypx^^. The identity matrix is denoted by I, 
Cfc is the unit vector such that {e-k)^. = 1 and {&k)i — 0, k I, and we also define vector 1 as (1)^ = 1, where in all 
three cases the matrix/vector dimension follows from the context. By {x,y) we denote the inner product of a;,y G C^. 
Without introducing ambiguity, we do not differ in the notation between random values and deterministic values. The mean 
of a random matrix (variable) X G C^^'^ is denoted as E{X). 



A. The parallel channels 

The set of parallel channels is denoted as A/" = {1, . . . , A^}. Let = {xki, ■ ■ ■ , xun) G 
be a random vector grouping the independent (zero-mean) symbols of user k E JC transmitted 
over the channels n E J\f equidistantly, at distance Tg. Then, the sampled signal of user k E JC 
received over the parallel channels can be written as = {yki, ■ ■ ■ , UkN), with 

where hk = {hki, . . . , hkN) G collects the path coefficients between the transmitter and 
receiver of user k E JC on channels n E M and = {uki, . . . ,nkN) ^ is a random 
vector which contains (zero-mean, independent of x^) Gaussian noise variables perceived at 
the receiver of user k E IC on channels n E M, where we assume = Ednkn]"^) > 0. The 
transmit power allocation to users and channels (in short, power allocation) can be written as 
P = (pi, . . .,pj^)' E M+''^, where vector = {pki, . . . .Pun) is such that pkn = E{\xkn\^) 
is the transmit power allocated to user k E IC on channel n E M. 

Let A = (tti, . . . , a^y E M;^^^ denote the sharing matrix of the channels among users 
such that flfc = (ofei, . . . , OfcAf) collects the relative fractions of time which are assigned to 
user k E K. for the exclusive access to channels n E M . Thus, as in practice the operation 
time is partitioned into frames of some fixed duration T ^ Tg, the collection of times Ta^ 
is reserved for user A; G /C for the exclusive access to the respective channels n G N within 
each frame. The set of allowed sharing matrices of the parallel channels takes the form 

^(r) = {AGM:^'^^: ||afc||i<rfc,A;G/C, a^,) < 1, n G AT}, r G R:^+. (1) 

k&K 

According to the first constraint in ([T]), a predefined vector r = (ri, . . . ,rx) G with 
||*"|li < N, is such that rk/N represents the fraction of the set of parallel channels which 
is assigned to user k E JC over time (over each frame). It proves useful in the remainder 
to introduce also R E R^""^ such that (i?)^,^ = 0, k ^ I, and (i?)^^ = r^, k,l E JC. For 
instance, under r = ^1 any user is assigned an equal 1 //^-fraction of the ensemble of 
parallel channels over time (over each frame). The second inequality in the definition ([T]) 
models then the obvious constraint that the aggregate time of exclusive uses of a single 
channel n E M hy the users k E JC does not exceed the total operation time (the total 
duration of each frame). Currently, the most celebrated instance of the considered parallel 
channels is the multi-tone/multi-carrier channel accessed by multiple users, as considered 
e.g. in [20], [21]. In this case, ak groups user's k E JC relative times of exclusive uses of 
carriers n E M and r^/iV represents the fraction of the multi-carrier spectrum which he is 
assigned over time [22]. 

Given a sharing matrix A E A{G, r) under use, we assume an arbitrary set V{A) of 
allowed power allocations, requiring merely that 

V{A) D {S,{0) n M;^^^), for some e > 0, e = e(A), (2) 

Such condition means, broadly, that all power allocations which are sufficiently small for 
the used sharing matrix are allowable. In particular, assuming frames of duration T, we can 
take either of the sets 

P(A) = {PGM^^^:5^r(a,,p,,)<i?}, A E Air), (3a) 

P(A) = {PGM^"^:T(afc,p,)<^fc, k E JC}, AEA{r), (3b) 

for some E, Ek > 0, k E JC, which mirror the limitations of energy per frame as a crucial 
constraint in current and future wireless communication systems [23], [24]. The latter set 
corresponds to conventional limitations of energy per frame in a multi-user Gaussian channel 
with user energy per frame budgets constrained hy E^, k E IC. The first set models the 



possibility of energy coordination among all users under the joint energy per frame budget 
constrained by E. This is the case, for instance, when the considered parallel channels are a 
means of representation of the orthogonalized broadcast channel which applies, suboptimally, 
single-user precoding instead of multi-user precoding [25], [26] (for the combination of 
parallel channels and the broadcast channel see also [27]). 

Complementarily to frame energy constraints it is sometimes desired to account for transmit 
power constraints at any time in a frame. In analogy to ([3]), under limitation of transmit 
power of any user k E JC hy Pk > and under the joint transmit power budget of all users 
constrained by P > we take, respectively, 

= P = {P G M^^^ : |||P|||i < P}, (4a) 
V{A) = v = {Pe M^"^ : \\p,\\, < Pfc, ke /C}. (4b) 

It is interesting to note that transmit power constraints at any time within a frame make 
the set of allowable power allocations independent of sharing matrix A E A{r) under use, 
which will be of key importance at several points in the remainder. 

B. The user performance 

For any user k E IC accessing the parallel channels, we assume a general vector-valued 
performance/ QoS function 

where we have fk{p) = (/fci(Pi), • • • , /fcjv(p7v)), with p ^ fkn{p), p>0,nEM. Function 
fkn expresses the performance of user k E IC on channel n E M, as a function of power 
allocated to channel n E M, when the user accesses this channel exclusively throughout the 
operation time. We restrict us to nonnegative QoS functions 

/fc(p)eM^^, pER'l, kElC, (5) 

and to avoid later technical queemess we assume that -^fkip) > 0, k E JC, for p E 
S'e(O) n and some e > (that is, performance functions are componentwise Frechet- 
differentiable and increasing at least for sufficiently small power allocations). 

Due to ([5]) and the assumed independent symbols of a user on each one of the parallel 
channels, it is reasonable to consider 

(a,p) ^ (ajkip)), (a,p) eR^ x R^, \\a\\, < r^, 

as the performance/QoS metric of user k E JC. Such metric represents the aggregate perfor- 
mance achieved by user k E IC on the entire channel ensemble, throughout the operation time 
(throughout each frame), as a function of powers allocated to channels n E M and relative 
time fractions of exclusive channel uses. We refer to a value of the performance metric of a 
user, for some A E A{r) and P E V{A), as user performance under policy (A, P). 

Let any predefined requirement/expectation of user k E JC with respect to the user per- 
formance be denoted as 7^ > 0, A; G IC. Then, min^gyc ^"'^'^^'^^''^'^ can be seen as the 
worst relative performance among the users accessing the parallel channels under a policy 
( A, P) G A{r) X V{A). Hereby, we implicitly assume that a smaller user performance 
implies a worse perceived service quality at the corresponding receiver. Such assumption 
complies with the nature of the very most QoS functions used in communications and 
information theory, but does not necessarily require strict componentwise increasingness 
of fk, k E JC. We give a few celebrated examples of such performance functions. 



Example 1 (Symbol decoding reliability): Let user k E IC access channel n E M and use 
uncoded constant-envelope modulation. Then, the achieved probability of error-free symbol 
decoding is 



p\h 



12 



with Q denoting the Marcum Q-function, M denoting the constellation size, and c > as 
some constant (e.g., c = 2 for binary Phase Shift Keying or c = 1 for binary Frequency Shift 
Keying) [28]. By ^ and the uniform symbol distance T^, the map (a^, p^.) ^{a^, fk{Pk))^ 
(A, P) E A{r) X V{A), expresses the aggregate (over channels n E M) average number of 
error-free decoded symbols of user k E K, in a frame as a function of policy. 

Example 2 (Mean square detection reliability): If the receiver of user k E K, utilizes the 
Minimum Mean Square Error (MSE) receiver and the user accesses channel n E H, then 
the achieved MSE can be expressed as — ^^^^ [29]. As a consequence, 

can be regarded as a kind of symbol detection reliability in the mean square sense. Thus, 
given (|7]), (ak,Pk) ^ ^{o-k, fk{Pk))' (A, P) E A{r) x V{A), describes the mean square 
detection reliability of user k E JC, aggregated over all symbols received in a frame, as a 
function of policy. 

Example 3 (Spectral efficiency): Let the modulation constellation size of user k E K, 
which accesses channel n G A/" be constrained by M. The spectral efficiency, in the sense 
of maximum number of reliably decodable bits/nats per symbols under given modulation 
constellation^, is not expressible analytically but can be approximated by a function 



r / \ rP\hkn 
fkn{p) = g[ n- 



2 



such that the map x t— > g{x), x > 0, is nondecreasing, and g{Q) = and gi^x) = log2 M, 
X > xo, for some a^o = xo{M) [8]. Then, it is easily seen that map (ak, Pk) ^ ^(^fc, fk{Pk))^ 
(A, P) E A{r) X V{A), describes the achievable number of reliably decoded bits/nats in a 
frame, as a function of policy. 

Example 4 (Capacity): If user k E KL utilizes the Maximum Likelihood (ML) receiver, 
then 

fkn{p) = log(l + (8) 
^kn 

represents the (information) capacity^, achievable by user k E JC when accessing the channel 
n E M, that is, the overall maximum number of reliably decodable bits/nats per symbol. Thus, 
given dSl), the function {ak,Pk) ^ ^(o-fe? /fc(Pfc))' (A, P) E A{r) x V{A), corresponds to 
the achievable (under Gaussian codebook) number of reliably decoded bits/nats per frame. 

With the given assumptions on user performance, the maximum attainable performance of 
the worst-case user accessing the parallel channels among policies from A{r) x V can be 
expressed as 

. {akJkiPk)) .o^ 
max mm . (9) 

iA,P)GA{r)xV{A) fce/C 7fc 

According to the common understanding of fairness in various multi-user channels, see e.g. 
[!]> [2], [11], [14], we refer to ^ as the max-min fair performance (in/of the considered par- 
allel channels), and we say that a pair (A, P) = argma,X(^A,p)eA(r)xViA) min^gAC ^"''"^J^^^''^'^ , 
is a max-min fair policy, which is, in general, not unique. 



'obviously, the notions of spectral efficiency and capacity are meaningful only when the duration of the considered 
channel access is sufficiently long, in the sense Tukn S> Ts. 



III. Graph of parallel channels sharing 

For any sharing matrix A E A{r) we define an undirected graph of parallel channels 
sharing, in short a sharing graph, which is induced by A. For the definition, recall that any 
graph is a pair, say G = {K,,8), where K, is the set of graph vertices, and 8 is the set of 
edges; any edge is represented by a pair {k, I) E S such that k,l E IC are the vertices which 
are joined/connected by this edge (are adjacent) [30]. 

Definition 1: For N > K and any A E A{r), a corresponding sharing graph G = G{A) 
is such that G = (/C, S) where {k, I) E £, k ^ I, if (a^, a/) > 0. 

The proposed induction of a sharing graph G by A is a version of orthogonal graph 
labeling from [31], which further differs slightly from the original concept of orthonormal 
representation of a graph in [17]. Precisely, a (not necessarily nonnegative) matrix A E M^^^ 
is referred to as an orthonormal representation of graph G = G{A) = (/C, £), which we write 
as A G A'^{G), if {ak,ak) = 1, k E JC, and {k,l) E £, k ^ I, whenever {ak,ai) ^ 0. By 
Definition [TJ any two vertices k,l E JC, k ^ I, of the sharing graph are adjacent if some of 
the parallel channels are shared by users k, I, where a channel is said to be shared by some 
two users if both users access this channel exclusively some fraction of time. The converse 
is also clear: If two nodes k,l E JC, k ^ I, of the sharing graph are nonadjacent, then no one 
of the channels n E M is shared by users k, I. 

It is readily seen that, for any fixed A E A{r), an induced sharing graph G = G{A) is 
in general not unique. Besides this, the graph characterization of parallel channels sharing 
provides merely the information on the topology, or combinatorics, of sharing relationships. 
Thus, given any graph G = {JC,£), different sharing matrices induce G as a sharing graph 
and we can group them in the set 

A{G, r) = {AE A{r) : G = G(A)}, r E M^+. 

The illustration is provided in the following example. 

Example 5: Consider parallel channels as a multi-tone/multi-carrier channel with N = 7 
tones accessed hy K = A users in the proportions r = (2,2,1,2). Let the sharing of the 
tones be described by the sharing matrix 

/ 0.2 0.5 0.5 

1 0.5 0.5 

0.3 0.3 0.3 

\ 0.5 0.7 0.7 

Then, the three possible sharing graphs G = G{A) are depicted in Fig. [H with the graph 
on the right hand side as the sharing graph with the minimal number of edges. On the other 
hand, for G as the minimum sharing graph from Fig. [H the set A{G, (2, 2, 1, 2)) of sharing 
matrices inducing it includes, in particular, column permutations of all matrices of the form 

ai2 ai4 \ 

024 ^25 
^32 ^33 ase 037 
042 ^43 046 ^47 / 

with ttij > 0, 1 < i < A, 1 < j < 7. 

A special role in our considerations of sharing graphs is played by the subgraphs called 
cycles. 

Definition 2 ([30], [32]): A cycle of length M in a graph G = {JC,£) is a sequence of 
distinct graph vertices ki E JC, 1 < i < M, which satisfy {ki, /ci+i) E S, 1 < i < M — 1 and 
{kM,ki) E £. 
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G = (/C,f) G = (/C,0 G={K:,£) 

Fig. 1. Three sharing graphs of the multi-carrier channel with A' = 4 users considered in Example 5; such graphs are 
induced, in the sense G = G{A), by the matrices A from Example 5. 




g = [k:,£) 

g = {k:,£) 

Fig. 2. Two exemplary sharing graphs of parallel channels with K = 7 and K = 12 users with marked exemplary cycles 
of length 3 (dotted edges), 4 (dash-dotted edges) and 5 (dashed edges). 

In simple words, a cycle represents a simple closed path in a graplJl. Note that the length 
of a cycle is the number of edges, or equivalently vertices, constituting the cycle. As an 
illustration, in Fig. [21 particular cycles are emphasized in two exemplary sharing graphs. A 
cycle of a sharing graph has an easy interpretation in terms of sharing policies: A cycle of 
length, say, M corresponds to a chain/sequence of M users accessing the parallel channels 
such that any pair of subsequent users shares some channel and the last user shares a channel 
with the first user. 

In Fig. [3] we show examples of M -partite sharing graphs for M = 2, 3, 4. As can be seen 
from the figure, such graphs contain only edges between some disjoint vertex subsets: The 
vertex set of an M-partite graph G = {}€,£) is divided into partitions /Cj, 1 < z < M, such 
that (k, I) ^ £ whenever k, I G /Ci, 1 < i < M. It is easily deduced that an Af -partite graph 
can not contain any cycle longer than M. An M-partite sharing graph is induced by sharing 
policies of parallel channels which distinguish M classes of users with the property that 
users within one class are not allowed, or not able, to share any channels over time. Such 
constraint is likely to be imposed by traffic processing regulations and/or the implementation 

"In the context of undirected graphs, some works prefer the notion of a circuit to the notion of a cycle used here. In 
such a convention, the cycle is understood as the analog to the circuit in directed graphs. 
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G = (/C,5) G = (/C,5) G = (/C,5) 

Fig. 3. Three exemplary M-partite sharing graphs of parallel channels with K = 5,7, 12 users and for M — 2,3,4, 
respectively. The graph on the left-hand side is a sharing graph of the multi-carrier channel from Example 6 and is induced 
(in the sense G — G{A)) in particular by the matrix A from Example 6. The graph in the middle can be a sharing graph 
of the multi-carrier channel considered in Example 7. 



effort, as is illustrated by the following examples. 

Example 6: Consider parallel channels as a multi-carrier channel with = 9 carriers and 
K = 5 users accessing the carrier set in the proportions r = (2, 2, 0.5, 1.5, 3) and having a 
common transmitter. Let the traffic of users 1,2 G /Ci be the real-time traffic, like voice or 
multimedia, while users 3, 4, 5 G /C2 transmit and receive so-called background traffic, such 
as file transfer, signaling or system information. From the viewpoint of percepted QoS and 
traffic processing complexity, it may be required to assign a carrier to real-time traffic for 
a large relative time fraction, say, no less than half of the total time, so that no carrier can 
be shared by two real-time users. Due to the processing effort, it may be also undesired to 
share a carrier between multiple users carrying the minor background traffic. These conditions 
enforce that the single carriers are either user-specific or carry mostly the real-time traffic 
of users JCi with some add-on background traffic of users /C2 from time to time. Thus, an 
exemplary sharing matrix can take the form 



G ^(3,2,0.25,1.25,2.5), 



which induces the bipartite sharing graph G = G{A) on the left hand side of Fig. [3l 

Example 7: Let a multi-carrier channel with = 13 carriers accessed by = 7 users 
in the proportions r = (3, 2, 2, 1, 1, 3, 1) be the considered parallel channels with a common 
receiver. Let the users be grouped in classes 1, 2 G /Ci, 3, 4, 5 G /C2 and 6, 7 G /C3 such that for 
any two users within one class the difference between their propagation times is larger than 
some critical propagation time difference (note that such classification is not always possible). 
Then, the sharing of a carrier between two users from one class can be undesired due to the 
required effort of time- and frequency synchronization to come up with the propagation time 
difference. This implies that the carriers are either user-specific or shared only across the 
classes /Ci, /C2, A^s, and that a particular sharing matrix A G ^(3, 2, 2, 1, 1, 3, 1) can induce 
the 3-partite sharing graph G = G{A) in the middle of Fig. [3l 
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As shown in the remainder, the description of the channel sharing topology by a sharing 
graph plays a key role in the problem of ensuring max-min fairness 



A. Selected algebraic graph characterizations 

We make use of the description of a graph by its so-called feasible matrix, which is a 
symmetric matrix indicating an edge by a nonzero entry [33], [34], [35], [36]: The set of 
feasible matrices of a graph G = {}C,S) is 

C{G) = {C= (cm) e §^ : cm ^ 0, A; ^ /, iff (A;, /) G S}. 

Given G, another set of interest here is parameterized by a vector v E of its vertex 
weights and can be written as 

B\G,v) = {B = %i) e §^ : hi = {vkViY^, {k,l) i£ork = I}. (10) 

For V = 1, this concept has its origin in the seminal work [17] where it was used in an 
approach to the problem of graph capacity. The generalization to the case v G was 
provided later in the study of relaxations of the vertex packing problem [18]). The graph 
capacity problem, as the key problem of zero-error information theory, remains still unsolved 
in the general case [19]. The set (flOl) is, however, a central element of the concept of so- 
called weighted Lovasz function, which in unweighted form (i.e. for v = 1) represents a 
general upper bound on graph capacity and is equal to the capacity for a certain graph 
class, including e.g. self-complementary graphs with vertex-transitive automorphism groups 
[17], [30]. Precisely, the weighted Lovasz function (later, simply Lovasz function) of a graph 
G = {IC,£) is defined as the map 

(G,v) ^e°(G,v) = min max , v eR^, (11) 

cGlR^:||c||2 = l 

with A'^{G) as the set of orthonormal representations of G, and it has the following property. 
Proposition 1 ([18], [31]): For any graph G = {IC,£), we have 

e\G,v)= min XmaxiB), v eR^^. 

In [18], [31], [37] one can find further interesting characterizations of the Lovasz function. 
A similar set which we make use of is 

B\G,v) = {B = {bM)eS'' -.hiyivkVi)-^, {kj)^£ ork = l}, (12) 

for any graph G = (/C, £) and v E M;^ is a vector of its vertex weights. For the case v = 1, 
the concept of the set (fT2l) is known from the considerations on the Delsarte bound, or 
Delsarte number, in [38], and the generalization to the case v E M;^ is straightforward. The 
unweighted (i.e. for v = 1) Delsarte bound was proposed within the algebraic framework 
of coding theory in [39], as an upper bound on the cardinality of an -clique, ^A C 
{1, . . . , A/}, in an association scheme with M associate classes denoted here as (/C, {£i}f£^). 
As the notion of association scheme is only loosely related with our topic, we refer here 
to Appendix |B] for an outline of the theory. One can relate the Delsarte bound for an 
association scheme (/C, {£i}fi-i^) to the graph G = {IC,£), Lifl-^£i, i.e. the graph whose edge 
set corresponds to the union of associate classes: The unweighted Delsarte number upper 
bounds the independence number of such graph [38] and the weighted Delsarte number (later, 
simply Delsarte number/bound), denoted as map 9^, has then the following characterization. 
Proposition 2 ([38]): For any graph G = {}C,£), we have 

e\G,v)= min XmaJB), v eR^. 



Further formulations of the Delsarte number which are direct analogs of the original char- 
acterizations of the Lovasz function from [17] can be found, partly without proof, in [40] 
and [41]. In particular, the authors apply the framework of graph Laplacians and identify 
the Delsarte number with the so-called ci-function of a graph and the Lovasz function with 
a related version of it. Similar characterizations of the Lovasz function and Delsarte bound 
and their properties in terms of edge orbits are studied in [42]. 

For our purposes, we define two further sets of the type (flOl ). ([T2l) and two related graph 
functions in the spirit of Propositions [T] and [2l First, we associate with a graph G = (/C, £) 
and a weight vector v G M;^ , the set 

B\G,v) = {B = (bki) e §^ : bki = (vkVi)-^, (kj) i£ork = l, 

hki < (vkVi)^, k,l e /C}. 
In analogy to Proposition [H we define for any graph G the map 

{G,v)^e\G,v)= min A„,,,-(S), v eR"^. (14) 

BgB2(G,ii) 

Second, also the set 

B''{G,v) = {B = ihui) e §^ : hi = {vuvif^, {k,l)^Sork = l ^^^^ 

vhvl' -B + \„a.{B)I e P^} 

associated with any G = {}C,£) and v E proves to be of key use in the remainder. By 
the definition of the class of completely positive matrices in M^^^ (Appendix lA]), the 
latter condition in (fTSl) can be written equivalently as 

v^v^' - B + X,„a,{B)I = V'V, for some VeMf^-^,MeN. (16) 

Furthermore, it is worth noting here that the condition v^v^ — B + \,„ax{B)I G P^^ is 
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implied by a slightly stronger requirement that v^V^ — B is included in the closure of P : 
This is an immediate consequence of \max{B)I = J2keic V ( ) ^fc V ■^max {B)e',^ and the 
characterization in Appendix [Al By analogy to Proposition [H for any G = {}C,£) we define 
a further map 

{G,v)^9'{G,v)= min A,„«,(B), v eR"^. (17) 

The relations between (flOl ). ([T2l) and the proposed sets ([T3l) . ([T5] ) are readily seen. It is 
immediate that v) C B^{G,v) and that the second condition in ([T3]) can be written 

as v^v^ —Be R^^'^ . Thus, by the form ([T6l) of the second condition in (flSl) and by 
inspection of (flOl) and ([T3l) . it can be seen that 

B\G,v) c B\G,v) c B°{G,v) c B\G,v), 

and thus 

0^{G, v) > e\G, v) > e%G, v) > e\G, v) (18) 

for any G = (/C, £) and v E R^ on hand. 

It is worth noting here that, given v = 1, the classes of matrices (flOl) . ([T2)) generalize the 
set of so-called (1, 5)-adjacency matrices of graph G introduced in [43]. Any (1, 5)-adjacency 
matrix is further affinely transformable to a Seidel adjacency matrix [32]. 



B. Some relations of the characterizations 

The algebraic graph descriptions introduced in Section IIII-AI have some simple properties 
which turn out to be central to our results. The first lemma below can be partially deduced 
from the proof of Theorem 3.5 in [35]. We give the proof for completeness and refer to 
Appendix |A] for the notions related to the set of completely positive matrices P^^ , such as 
the cp-rank. 

Lemma 1: Given any graph G = {}C,S), we have 



and the cp-rank satisfies 



<P{C) < ^^^^^\ CGC(G)nP^. 

Proof: Associate any edge (A;, I) E S with e = e(/c, /), 1 < e < and let B = (bki) G 
rJ'''^' be defined a^ 

Kk>0, Ki>0 iff e = e{kj). 

Then it is readily seen that {BB')^^ > 0, A; 7^ / iff {k, I) e S, so that C = BB' satisfies 
C e C{G) and, by Definition [6l we also have C G P-^^. This proves C{G) n P-^^ ^ 0. 

According to the known bound on cp-rank, see e.g. Section 1 in [33], if additionally 
C e P^^, then we can find B e such that C = BB' for some N < K{K + l)/2, 

which completes the proof. ■ 

The lemma says essentially that the set of feasible matrices includes a completely positive 
matrix for any graph on hand and any such matrix remains completely positive if all diagonal 
elements are replaced by the largest eigenvalue. Furthermore, for any graph with K vertices, 
e.g. a sharing graph of parallel channels accessed by K users, any of its completely positive 
feasible matrices has a cp-rank no larger than K{K + \) /2. The latter bound on the cp-rank 
is the best known, but likely not the best possible bound [33]. 

Lemma 2: Given any graph G = (/C, S) and v G M^, consider the set 

Vi{G, v) = {v^2v^^' -B + X„u,,{B) : B G B'{G, v)}, i = 0, 2, 3. 
Then, we have 

Vo{G, v) c Ug'cgC{G'), V2{G, v) c UG"cGC(G")nM^^^, V,{G, v) c UG'cGC(G")nP^, 

where G' C G denotes that G' = (/C', S') is a subgraph of G in the sense that /C' C /C and 
S' C S. Moreover, we have 

V2{G, v) = Vs{G, v), equivalent to B'^{G, v) = B^G, v), 

if either K < 4 or G has no odd cycles longer than 4. 
Proof: Let any C = {cm) G R^""^ such that 

C = v-^v'2' -B + X,„„,{B), (19) 

for an arbitrary B G B'{G, v), i = 0,2, 3, be given. Then, by the definitions ([JO]), ([H]), ([T5l) 
we have 

Cki = 0, k^l, if {k,l)iS, 
but also Cki = 0, A; 7^ /, if [k, I) G JF, for some T '^E, where JF = is allowed. This implies 

Cfci^O, ki^l, iff {k,l)^E', 

'in the particular case h^h = = 1 iff e = e{k,l), matrix B represents the so-called incidence matrix of graph G 

[30], [32]. 



with E' = E\T . Thus, for any matrix C we have C G C(G") for some subgraph G' C G, 
with G' = {IC',£'). 

If now i = 2, then it is evident by the definition ([H]) and b>0 \max{B) > that C eB^ 
(see Definition [5]). By the result in [44] this implies C G whenever A' < 4 (see Appendix 
lAl) . Since C G C(G') for some G" C G is proven for any (fT9l) . we have further by Theorem 
3.1 in [34], or by [45], that C G holds also if G has no odd cycles longer than 4. By 
the definition ([T5] ). this completes the proof. ■ 

An implication of the lemma is that any matrix C G Vi{G,v), v G M;^, is a feasible 
matrix of some subgraph of G. Further, any matrix C G V2{G,v) is completely positive 
whenever either a graph G with no more than 4 vertices is considered or when the maximum 
odd cycle length in the graph is no longer than 4 edges (the existence of some completely 
positive matrix C G 1^2 (G, v) is ensured already by LemmaH]). In particular, any such matrix 
is completely positive for G as a sharing graph if the parallel channels are accessed by no 
more than 4 users, or if there are M < 4 classes of parallel channels users, where channel 
sharing within a class is not allowed/possible due to restrictions on implementation or QoS. 
Recall that such parallel channels are illustrated by Examples |6l |7] and their graphs are given 
in Fig. |3l 



IV. Upper bounds on max-min fair performance 

In this section we derive several upper bounds on the worst-case user performance in the 
considered parallel channels. According to our performance model, an upper bound represents 
an optimistic case, i.e. a better value of user performance than the upper bounded one. The 
bounds in this section are not proven to be tight and thus, are not very interesting when 
considered alone. They become, however, interesting and lead to the central conclusions of 
this work when considered together with the lower bounds from Section |Vl 



A. Upper bounds 

In the following Proposition, a policy-specific bound on the worst performance within the 
user population is proposed. 

Proposition 3: Given N > K, any G = {)C, S) and any {A, P) G A{G, r) x V{A), 
r G M;^, we have 

inin 5 < — ' ^, r , « = 0,1,2, (20) 



with w G such that 



wk = ^, keJC, (21) 



and where we defined 

fiA,P) = {feR^:{akJ)>{akJk{Pk)),ke}C, for some AG^(G,r)}. 

Proof: Given any A G A{G, r), f eR^ and any P G P(A) such that {ak, fk{Pk)) ^ 
0, k E K. (which by our assumptions in Section HI] exists) , let us define Z = (zi, . . . , zk)' G 
R-^^^^ such that 

Zk = Zk[ak) = Jj:^-^f - —f — fc G /C, 

V(/,/) {akJkiPk)) 

''This is readily seen by tlie feature tliat for any C — BB' we fiave \mm{C) — niax^gjjK.||3.|| x'B{x'B)' > 
e',B{e',BY = {BB'),, > 0, k e IC. 



with an arbitrary A = (ai, . . . , aj^)' G R . Then, we have 

{Zk,zi} = y/WkWi - — — — +- — — TT^' ^'^^''C, 

{o-kJkiPk)) {aiJiiPi)) {akJk{Pk)){aiJi{Pi)) 

(22) 

for any A E M^^^. Let now / G JF(A,P), and note that then we can find a particular 
A G A{G, r) which satisfies 

{akJ) = {akJk{Pk)), keIC (23) 

(in fact, the system (l23l) has always a solution A G R-^^^ and since by / G ^{A, P) there 
exists some A G A{G, r) such that {cik, f) > {ak, fk{Pk))^ k E JC, it is implied that A < A, 
and thus A G A{G, r)). When A G A{G, r) satisfying ([23]) is taken in we yield 

{zk, zi) > -y/wkWi, kJeJC, (24) 

where in particular (zfc, Zk) = —Wk + "'(i" /"(Pfef)'^"^^ ' ^^'^'^^ Definition [T] (/c, I) ^ £ 
implies (a/,., ai) = 0, also 

{zk,zi) = -y/w^i, {kj)^£, k^l. (25) 

Thus, by the definition (fT3l) . we can write 

- B = ZZ' - (/, f)G{A), for some B G B'^{G, w), (26) 

map A ^ G(A), A G M'^''^ 

r'iiakJkiPk))'^ 



where the map A t-^ G(A), A G , follows by the definition of w as (G(A))^^ 



A; G /C, and (G(A))^^ = for k,l e JC, k ^ I. Feature ([261) implies then 

max -By ZZ', (27) 

fce/c {akJk{Pk)? 

since by such definition of G and by the property ^"^^'^''^ = < 1 we have (G)^^ < 



2 

maxfcgyic 7 — Yi — TT^' k E IC. For the particular S G B'^iG, w) in (|26l) we have then 

ijk \Pk)) 



max f >%^. /e^(AP) 
fcex {akJk{Pk)r if J) 



so that for i = 2 the result follows by the definition (fT3l) . For the cases i = 0, 1 the proposition 
follows then from the definitions (flOl ), (fT2l) and from the property (fTSl ), which completes the 
proof. ■ 

By the proposition, the worst squared user performance achieved under any policy (A, P) G 
^(r) X P(A) in parallel channels is no better than the ratio of the least 2-norm achieved 
among the vectors within the set JF( A, P) and the function 6*^ evaluated for a sharing graph 
induced by A and for the vector w such that (|2TI) . According to (fTSl) . when such value of 
9"^ is replaced by the Lovasz function value or Delsarte bound value assumed by the sharing 
graph and the vector w, the bound from Proposition [3] is loosened. Since 7^ is a predefined 
performance requirement and the fraction of the channel set J\f assigned to user k E JC 
over time, w can be interpreted as the vector of squared user performance requirements 
normalized by assigned channel fractions. 

It is readily seen that JF( A, P) is the set of values of performance function^ which 

• are equal for any user accessing the parallel channels and, 

• for some sharing matrix which induces the same sharing graph as A (i.e. under fixed 
sharing graph), attain user performance no worse than under policy (A, P). 

^The value of the performance function /fe, fc £ /C, is a vector in M.^ and shall not be confused with the user performance. 



see our performance model in Section HI-BI 



Thus, in some sense, J^{A, P) can be seen as a set of dominating values of QoS functions 
for the policy {A, P). Note that a QoS function value / e ^{A, P) may be not achievable 
by an allowable power allocation from Vi^A), as such value leads to a superior multi- 
user performance under the penalty of being equal for all users. The set P) is not 
a polyhedron for a general {A, P) E A{r) x V{A). Nevertheless, for any given policy 
(A, P), T{A, P) contains the polyhedron 

{feWl: {ak, f) > {ak, fkM), k e JC} (28) 

and its further polyhedral subset {/ G M.^ : / > fk{Pk)^ k E IC} which depends merely 
on P. Both polyhedra give rise to obvious simplifications of (|20|) : In particular, for any 
(A,P) E A{r) X V{A) and for U^^{P) = {maXk(.K: fki{Pk), ■ ■ ■ ,^^^k&K fkN{Pk))' we 
have 

. {akJk{Pk)? ^ (/max(f ),/max(f )) {Y.keK fk{Pk) .Y^keK fk{Pk)) ■ _ n 1 9 

k^K 72 - 6i{G{A),w) - 9i{G{A),w) ' ' ' ' 

with w such that (|2T1) . Thus, given any policy in parallel channels, the worst squared 
user performance can be no better than the squared 2-norm of the channel-wise maximum, 
respectively sum, of performance functions of users divided by the function 0"^ (or the Lovasz 
function or the Delsarte bound) evaluated for the induced sharing graph and the vector of 
squared user performance requirements per assigned channel fraction. 

The technicality of the bound (|20|) lies in the structure of the optimization domain JF( A, P), 
while the weight vector w is easily interpretable. As Corollary [8] in Appendix O we prove an 
alternative version of Proposition [3] which simplifies the optimization domain in the bound at 
the expense of a more complex weight vector structure. The bounds from Proposition [3] and 
Corollary [8] yield the following implication on the max-min fair performance under given 
sharing topology of parallel channels. 

Corollary 1: Given N > K, any G = {JC,£) and r E IR+4., we have 

{akJk{Pk)? < "^i^/g^(G,r)( /, /) 

{A,P)(^A{G,r)xV{A)Tiic 7^ " 9''{G,W 

where w is such that (|27]) and where 



max, _,.,min ^ < "^Cr ^ V , ^ = 0,1,2, (29) 



:F{G,r) = {fER'l:{ak,f)>{akJk{Pk)), k E IC, for some AEA{G,r)}, 
with 

/4 ox . {o-kJkiPk)) 

(A,P) = arg max mm k . 

{A,P)eA{G,r)xV{A) keIC 7^ 

By Proposition [3] it is evident that J^{G, r) is equivalent to the set of dominating QoS function 
values P), where (A, P) E A{G, r) x V{A) is a max-min fair policy under a fixed 

sharing graph G. As JF(A, P) contains the polyhedron (l28l) for P = P, we get the following 
loosened version of (|29l) . 

Corollary 2: Given N > K, any G = (/C, S) and r E vve have 

n,a. ^, J-^Jf^))' ^^ '-'-ff>o.r^(f-f\ ,^0.1,2, (30) 

{A,P)eA{G,r)xV{A) keK. 7^ 0\G,w) 

where w is such that (|27]) and where, with A defined as in Corollary [7] 

T{G,r) = {fER^:\\f\\,>\\fk{Pk)h, kEiC, PeV{A)}. 
Proof: First notice that for the policy (A, P) defined in Corollary \T\ we necessarily 
have llafell-,^ = r^, k E IC. Further, as for any A E A{G,r) such that = r^, k E IC, 



it follows that J^keic^k = 1, we can write the condition (/, 1) > (/^(Pfc), 1), k E JC, 
P G V{A), specifically as 

V(/,afe) >max(/z(p;), Vafc) = V(max/z(pJ,afc), AeA{G,r), PeV{A). 

fce/c keK keK. 

This further implies for a particular P = P that 

$^(afe,/) >$^(afc,/fc(p,)), AG^(G,r). (31) 

fce/c fee/c 

Let now A be defined as = afk{Pk) o ak o f ~^, k E IC, where /"^ = (l//i, . . . , V/at) 
and a > is chosen to ensure Wa^Wi < r^, k E /C, and thus A E A{G, r) (it is evident that 
any sufficiently small a satisfies such condition). For this particular A we have {a^, f) = 
a{ak, fk{Pk))^ k E }C, so that together with (|3T1) it is implied that a > I and on the other 
hand 

(ftfc,/) > (akJkiPk))^ kElC. 

Consequently, r) C r) which, by Corollary [H completes the proof. ■ 
The set r) includes all QoS function values, equal for all users, which are in the sum 
over all channels superior to any QoS function value achieved by an allowable (for some 
A E A{G, r)) power allocation. Thus, jF'{G, r) can be seen as a hull of any user dimension 
of the feasible QoS/performance set of parallel channels, which we define in analogy to the 
theory for channels with interference as [46] 

{(/i(Pi)),...,/i.(Pi.)):PeP(A)}, (32) 

(equivalently, yikeK.^{G,r) is a hull of the feasible QoS set). 

Corollary [H implies that a squared max-min fair performance under the condition of a fixed 
sharing graph G in parallel channels can never exceed the ratio of the minimum squared 2- 
norm within the hull !F{G, r) of any user dimension of (l32l) and the value of the function 6"^ 
(or the Lovasz function, or the Delsarte number) assumed by G and the vector w satisfying 

m. 

Consider now constraints on transmit power at any time (in a frame), as expressed e.g. 
by (l4al) . in which case we have V{A) = V, A E A{r) (allowable power allocations are 
independent of sharing matrices and sharing graphs). In such case it is readily seen that also 
set ^{G, r) is independent of the sharing graph on hand, i.e. !F{G, r) = P{r) regardless of 
G, and thus the bound (|30l) assumes a specific separated structure. Precisely, the max-min fair 
performance under a fixed sharing graph is upper-bounded by a ratio of a value dependent 
solely on this graph and a vector norm determined completely by the the attainable power 
allocations. Thus, (l30l) provides a separation between the influence of the combinatorial 
topology induced by the channel sharing policy via Definition [U and the impact of (the 
structure of) the set of allowable power allocations. The optimistic bound (l30l) . although 
looser than the one from Corollary [H proves in the next section to be particularly insightful, 
since a complementary pessimistic bound of the same type can be given. Again, recall that 
according to (fTSl ), Corollary [T] and (|30l) provide the tightest bounds when the extension 6"^ 
of the Lovasz function and the Delsarte bound is incorporated. 

Obviously, we can reformulate Corollary [T] and (l30l) for the max-min fair performance 
nonrestricted in term of the sharing graph. Precisely, 

(,A,p)eA(,r)xV{A) keic -fl 9'{G,w) e'{G,w) 



with G as the max-min fair sharing graph in the sense that A G A{G, r) (equivalently, 
G = G{A)), where now 

(A, i-')=arg max mm ^ (33) 

iA,P)eAir)xV{A) fce/c 7^ 

is the (graph-nonrestricted) max-min fair policy of the parallel channels. 

B. Relations to coding and zero-error capacity 

Relations of max-min fair performance in parallel channels to coding and zero-error 
information theory results are obtained in the setting 

— = 1, kelC. (34) 

Tk 

This can be assumed for a homogeneous user population, that is, if an equal fraction of 
the parallel channels is to be assigned (over time) to any user and all users have equal 
performance requirements. By the celebrated result in [17], the Lovasz function of G,w 
satisfies in such case 

e^{G,w) > e(G), 

where Q{G) = lim„_^oo ^/ct{G''^) represents the (zero-error) capacity of G; a expresses 
hereby the independence number of a graph and C?" denotes an n-fold concatenation, or 
power, of graph G [30]. The capacity interpretation of Q{G) originates from the fact that 
represents the maximum number of n-letter messages which will not be confounded 
when k E }C correspond to alphabet letters and any edge {k,l) E S models the (danger of) 
confusion of letters k, I [17]. As a consequence of Corollary [H (|30l) and the result of Lovasz 
we yield for any sharing graph G that 

2 ^ min_^g^(G,i)(/,/) min_^g^(^^^)(/,/) 



niax mm{akJkiPk)) < 7^77T^ < 



{A,p)eA{G,i)xV{A) keic' - Q{G) - e{G) 

In words, under a homogeneous user population accessing the parallel channels and under 
sharing graph fixed to G, the max-min fair performance never exceeds the minimum 2 -norm 
within the set ^{G, 1) (respectively, within the hull T{G, r)) divided by the square root of 
the sharing graph capacity. This means also that the max-min fair performance scales at most 
with the capacity of the corresponding sharing graph G = {)€,£), i.e. with the effective size 
of the alphabet needed for error-free communication of the letters /C where the letter pairs 
£ are confusable [17]. 

Given (|34|) . we have also the central relation of the Delsarte bound of G and the graph's 
independence number according to [38] 

e\G,w) > a{G), 

(recall that by ^ we have additionally 6°{G, w) > 9'^{G, w)). Thus, 



2 / min^g^(G,i)(/,/) min_fg^(c,^)(/,/) 



max mm{ak, fk{Pk)) < 77T^ < 



{A,p)eA{G,i)xV{A) kejc' a{G) a{G) 

which, with definition of the independence number, means that the ratio of min^gjr(G,i) {f, f) 
(respectively, minj,^.p(^Q ^-^{f , /)) and the maximum cardinality of a vertex subset of a shar- 
ing graph G such that no two vertices in it are adjacent upper bounds the max-min fair 
performance under fixed sharing graph. This implies that the max-min fair performance in 
parallel channels scales at most with the independence number of the sharing graph. 

We close the discussion of the upper bounds by pointing out two crucial issues. First, the 
given upper bounds on the max-min fair performance apply to the case N > K, i.e. to the 
parallel channel instances with the channel ensemble no smaller than the user population 



accessing them. Thus, the bounds apply to, in some sense, non-overloaded parallel channels, 
which allow the possibility of permanent (i.e. in each frame) access to a channel for any user. 
Second, the generality of the upper bounds has to be underlined. The bounds apply to any 
performance function for which the formulation of the max-min fair performance according 
to ^ is meaningful, that is, when a larger user performance implies a better perceived service 
quality level at the user receiver (Examples [HH]). 



V. Lower bounds on max-min fair performance 

The lower bounds on max-min fair performance presented in this section correspond to 
pessimistic values, in the sense that the max-min fair performance is guaranteed to be no 
worse. These bounds are analogs, or complements, of the optimistic bounds from Section 
HVl and together embrace the max-min fair performance in parallel channels. 



A. Some notes on matrix scalings 

The proposed bounds make use of some novel elements of the theory of matrix similarity 
and matrix scaling which are outlined in the following. Let us define a scaling of a nonneg- 
ative matrix by straightforwardly extending the idea of scaling of a square positive matrix 
from [47]. 

Definition 3: A matrix A E R;^^^ is said to be (r, c)-scalable, where r E and c E 
ra)iv ^ ij- ii^ii^ ^ ii^ii^ if there exist X = diag{X) E M+''-^ and Y = diag{Y) E R^""^ 



such that 



XAYl = r, I'XAY = c! . (35) 



The pair (X, Y) is then referred to as an (r, c)-scaling of A. 

Thus, an (r, c)-scaling of a nonnegative matrix collects scaling factors of rows and columns, 
in the form of two diagonal matrices, such that row sums grouped in r and column sums 
grouped in c are obtained under row-wise and column-wise scaling. A related notion which 
proves useful in later considerations is the set 

X{A, r, c) = {x = XI, y = Yl: {X, Y) is (r, c)-scaling of A E M^^^, (f, c) < (r, c)}. 

In words, X{A, r, c) consists of vector pairs which collect diagonal entries of those (f, c)- 
scalings of A e M;^''^ which are no larger tharH (r,c). 

Given predefined r E and c E it is obvious that matrices which are not 

(r, c)-scalable exist in M.^^ . Nevertheless, for any nonnegative matrix we can always find 
a scaling which leads to row and column sums no larger than the predefined ones. 

Lemma 3: Given any A E M.^^^ and any f E c E there exist r < r and 

c < c such that A is (r, c)-scalable. 

Proof: Let A = (ai, . . . , a^-)', with a/c E M^, k E JC, and define A = XA, where 
X = diag{X) is such that XI = x and 

Xk = T"^' kE IC. 
{cik, 1) 

Then, letting A = (ai, . . . , aj^)', we have Al = f, so that if I'A < c, the proof is 
completed. Otherwise, let A = AY, where Y = diag{Y) is such that Yl = y, with 

Vn = min , nEAf. 

Then, it is evident that 1' A < c! . Further, as [I'A)^^ > c„ for some n E M (by assumption), 
we have < y < 1, which implies also Al < Al = f and completes the proof. ■ 



Here and hereafter we refer to a scaling as larger/smaller than an other scahng if the obtained column and row sums 
are componentwise larger/smaller. 



The original characterization of a scaling (of a square positive matrix) was given in [47] in 
terms of a nonlinear program. The currently known descriptions of scalings of nonnegative 
matrices are mostly in terms of optimization problems, see e.g. [48] and references therein. 
In the following we provide a novel (to the best of our knowledge) characterization which 
extends the concept from [49]. 

Lemma 4: Let A G M;^^^ be {r, c) -scalable for some r G c G Then, if we 

define f = [r' 0)' G and if y E , satisfies 



V(^(t/) < (36) 



for the function 



^(^) = - J]log— 



z G 



and if X E is such that 

Xk = Xk[y) = . . . , /c G /C, (37) 

then {X, Y) such that XI = x and Yl = y is an (r, cj-scaling of A. Moreover, is 
satisfied if and only if y is a global minimizer 

y = arg min - V log (38) 

Proof: By the definition, we can write {Vip{z))^ = J^keK'^knTpH —, n E M, for 

any z E so that with (l37l) we have in particular for z = y that 

(Vv.(?/))„ = (A'x)^ - ^ = (A'X1)„ - ^, uEAf. 

yn yn 

This implies together with (|36l) that 

I'XAY < c. (39) 

Further, we have 

{XAYl), = Xk{Ay), = n, kEiC, 

by the definition dTT] ). and thus I'XAYl = I'c, since I'c = I'r holds by assumption 
(Definition |3]). Consequently, (|39| ) is satisfied only if XAY = c, and thus V^p^y) < only 
if Vip{y) = 0. To prove that the latter condition is equivalent to (l38l) . apply the transform 
V = log z, z E and then rewrite if with the properties of the logarithm as 



As r G and the map v ^ , f^s , v G M^, is known to be log-concave (see, e.g., [46], 

Chapter 6), it is immediate that v ^ '^{e") is convex for v E M . Thus, Vip{e'^) = is 
equivalent io w = argmin^gigiv '^{e"), which by the one-to-one setting w = logy gives (l38l) 
and completes the proof. ■ 
It is worth mentioning that function Lp from the lemma is multiplicatively homogeneous 
in the sense that 7^(2) = ^{az) for any z E M;^^ and a > (so that any minimizer (l38l) 
scaled by some a > is a minimizer of Lf as well). This is readily seen from the exponential 
transformation 



e 



ip(z) ^ nfc=i i^^)k 

1 Ln=l 



as used originally in [49], and from the condition ||r||^ = ||c||^. Furthermore, there is a 
surprising relation of function Lp to the throughput optimization under interference. Let 



us interpret z G as a transmit power vector of the user population M accessing the 
interference channel which has (A' 0)' G M^^^ as its interference matrix, defined in the 
usual way as e.g. in [46], [50] (this implies that the channel gains of N ~K users are zero). 
Then, by defining the Signal-to-Interference functions of users in the interference channel as 
z ^ SIR„(z) = j^^, n [50], we can write 

K K N 

(p{z) = - ^ rfc log SIRfc(z) - ^(cfc - Tk) log 2fc - ^ Ck log Zk- (40) 

k=l k=l k=K+l 

By this form, —Lp can be recognized as the weighted throughput function of the described 
interference channel with additional cost functions. When Ck > rk > 0, k E K,, such cost 
functions penalize logarithmically the excessive use of transmit power by the users. By 
the proof of Lemma HI the weighted throughput function (|40l ) is known to be convex as a 
function of the logarithmic power vector v = log 2, z G M^.,. (e.g. power allocation in dB). 
Lemma|4]and the above interpretation lead to the conclusion that (X, Y), with Yl = y and 
XI = X, is an (r, c)-scaling of an ((r, c)-scalable) A if y represents a power allocation 
which globally minimizes the penalized weighted throughput function (|4Ql) in the described 
interference channel and x is determined by y via (ITT] ). 
Finally, we need the following scaling-related function. 

Definition 4: Given r G let the map V /i(V), V G M+^^, be defined 

u(V) = max min {xkllnf' ■ 

{x,y)(iX{V,r,l) {n,k)€AfxK 

Such function represents the minimum squared geometric mean of pairs of diagonal entries 
of an (r, c)-scaling of a given matrix, achievable among all (r, c)-scalings no larger than 
(r, c). In the spirit of [51], we can regard /i as a (kind of) metric, or measure, of the entire 
class of such scalings of a given matrix. 

As the simple property of later interest, we observe that if the row and column sum vectors 
of V do not exceed (r, 1), i.e. VI < r and I'V < 1', then fi(y) > 1. 



B. Lower bounds 

Using Definition |4] we can formulate the following lower bound on the max-min fair 
performance in parallel channels under fixed sharing graph. 
Proposition 4: Given any G = (/C, S) and r G we have 

. {dk, fk{Pk)? ^ Ai(^T^) max^g^(^ ) (/, /) 
max mm ^ > max ; — r . 

(A,P)G.A(G,r)xP(A) fce/C 7^ BeBHC'w), VgM^x^: \nax[B) 

(41) 



VV' = X,„a:,'^{B){whwh' ~B) + I 



where w is such that d27]) . where 



N > (42) 



and where, with A defined as in Corollary [7] 

:F{G,r) = {feR'^:f = fk{Pk), keJC, for some PeV{A), 

. . {akj)\ 

t = arg max max mm „ - \ . 

Moreover, for a particular 

B = arg _ min \,„ax{B) (43) 



'We omit here the indication of the dependence on r, since it does not introduce any ambiguities in the remainder. 



(HIP implies further 

max mm 7. > max — ^- — , 

{A,P)&A{G,r)xV{A) keK. Yk VsR^^^: e^\G,W) 

VV'=X,„a.''^(B){whwi' -B)+I 

(44) 

Proof: Let any G = (/C, £") and any B E B^{G, w) be given, and let V = (vi, . . . , v^)' 
satisfy 

= X^ax^\B){w-^w-^' -B)+I, Ve M^^^, (45) 

where by the definition (fTSi) it is known that such V exists whenever = N(B) satisfies 
N > (l){\„,a.x~\B){w^w^' - B) + I). Defining = diag{W^) as (VT^)^^ = (lo^)^, 
k E /C, the right-hand side of (l45l) can be rewritten due to — S + \,„ax{B)I ^ as 

Xnu,x~\B){w^w^2' -B) + I = \,„ax~\B)W^^CC'W^^ + XX', (46) 

with any C = (c, . . . , c)' e R-^^^ such that (c, c) = 1 and with any X = {xi, . . . , xk)' E 
rKxn which satisfies XX' = -B + A„„„(B)/ and (ajfc)^ = 0, n > K, k E IC. Letting 
now > maxjif + 1, (l){\,„ax~'^{B){w^w^ — B) + /)}, we can find for any such X some 
C = C{X) E M^^^ satisfying 

(c, Xk) = 0, kElC, (47) 

so that (|46l) can be rewritten as 

X,„„x~\B){w^w"^' - B) + / = {\,„ax^HB)w"^C ± X)(A,„„.r^(S)i^^C ± X)'. (48) 



By (l45l) . it follows now that any vector tuple x^ + ^/ a "'fg) '^' k E IC, yielding (l48l) satisfies 



i.e., vector tuple + y y^^^c, k E IC, has the same lengths and mutual angles as any 

vector tuple v^, k E IC, yielding (l45l) . As a consequence, for any tuple v^, k E IC, satisfying 
(|45l) and for any c and cc^, k E IC, from (l48l) . there exists a rotation matrix Q E M^^^ (a 
real- valued orthogonal matrix with unit determinant) for which [52] 



Vk = Q{xk + Jy~^^^c), kElC. 

By orthogonality of Q we have {Qc,Qc) = {c,c) and {Qxk,Qxi) = {xk,Xi), k,l E JC, 
and (l47l) implies (Qc, Qtc/,) = 0, G /C. Thus, it follows now that any factor in (|45|) can 
be written as 

V = Jx,nax-\B)W-^C ± X (49) 



for some C = (c, . . . , c)' G R^""^, {c,c) = 1, and for some X E R^'""^ satisfying (gV]) 
(where X is such that XX' = —B + \,„ax{B)I). This further yields that 



and, by the Definition (fT5] ). also 

(-j^fc, = 1, A;G/C, (t;,, t;;) = 0, A; ^ /, if (A;,/)^£. (51) 



G /C, (52) 



By dSOj) we have 

\nax{.B) Wk 



(/,/) («.,/)^' 



for any / G chosen to satisfy , = c for the particular vector c in (l50l) . By Lemma 

V (/./> 

Hit is further implied that there exist (Z, Y) e M+'''^ x M^''^ which represent an (r, c)- 
scaling of iiV such that (r,c) < (r, 1): By setting Zl = z, Yl = y this means that 
we can take any {z,y) E X{RV,r,l), so that by (ISTI) and Definition [T] it follows that 
U = ZRVY satisfies U E A{G,r). Furthermore, writing U = (ui, . . . , ua')', we have 
then by the definition of w that 

Tfc 7l4 ^ ll min(fc,„)g/cx^(^fcyn)^ 

(53) 

which implies with (|52|) and Definition |4] finally that 

A„..(B) (54) 



f,{RV){fJ) - {u,JY 

Note now that (l54l) holds for any / E with J- — = c for the particular c in (l50l) and, 

V (/>/) 

by the assumption (|2l) and the assumptions with respect to fk. A; G /C, we can always find 
a particular / such that additionally / = fk{Pk)^ k E IC, for an arbitrary A E A{G, r) and 
for some P G V{A) i Consequently, it is further implied that 

— TTTTTr > mm m.&y.-^ — = max -^^^ — , fEJ^[G,r), 

fi{RV) - AeAiG,r)jeR^: fce/c {ak, fce/c (a^, /)2 

f=fk{Pk)MK., for some PgP(A) 

(55) 

where A is defined as in Corollary [T] and A = arg minyie^(G'^^) maxkfzjc "^"^p"- Thus, finally 
it is immediate that 

^""^^^ > min max- r^, /G^(G,r), (56) 



fi{RV){f,f) {A,p)eA{G,r)xV{A) keK {akJk{Pk)y 



for any B E B'^{G, w), for any V satisfying (|45l) and > max{_ft'+l, </)(A,„„;i~^(S)(if;2iu2 — 
B) + I)}. According to Lemmas [H [2l the latter condition is satisfied regardless of -B G 
B^{G, w) if N > K{K + l)/2. As is satisfied in particular for B such that A,„,„(B) = 
in.axBeB-^{G,w) \nax{B) = 9^{G, w), the proof is completed. ■ 
The proposition says that the squared max-min fair performance achieved in parallel 
channels under fixed sharing topology is guaranteed to be no worse than the maximum 
ratio of some two expressions. The denominator expression is the maximum eigenvalue of a 
matrix B from w), which is determined by the given sharing graph G and the vector 

of squared user performance requirements normalized by assigned channel fractions. The 
numerator corresponds to the squared 2-norm of a vector from the set ^{G, r) multiplied by 
the value of the metric ji of the class of (r, c)-scalings no larger than (r, 1), of a nonnegative 
factor of ^ 

R{\nax-\B){w-2w"^ -B)+ I)R. (57) 

By Lemma [2l the matrix (l57l) represents a particular feasible matrix of a subgraph of the 
sharing graph. Obviously, the looser bound (|4T1) is obtained by replacing the maximization of 
lji{RV) / \,nax{B), conducted over B E B^{G, w) and the factors of (l57l) . by the minimization 

'^Equivalently, by these assumptions, {/ £ : / = fk(pi^),k € IC, P G V{A)} has a nonempty intersection with 

the ray {/ G R+ : -yi — = c} for any c £ K^, and A G A{G, r). 
V {/•/) 



of the eigenvalue only. By Corollary \T\ and Proposition |4] we have now 

min^g^(G,^) (/, /) min_^g^(G,r) (/> /) 

— max 



\2 



e'^{G,w) BeB^Cw) XmaxiB) 

> max mm ^ > 

{A,P)(iA{G,r)y.V{A) keK 7^ (58) 

fi{RV) max f J f,f) ^{RV)maXj.^^,^ Jf, f) 

jy^g^^ Z. i J> i i 

BgB3(G,-lu), VGMf >^max{B) ~ 6^(0, w) 

VV'=X,„ar^{B){whw^ ~B)+I 

where RV in the outer lower bound denotes any nonnegative factor of (l57l) for the particular 
matrix (l43l) . achieving the value of the 6^ function (see (flTl)). While ^{G, r) was shown to 
be the set of dominating performance function values for some policy (A, P), set P{G,r) 
includes precisely those QoS function values which 

• are equal for any user accessing the parallel channels, 

• are attainable by some allowable power allocation (under some A E A{G, r)) and, 

• optimize the worst user performance under fixed sharing graph G and under QoS 
function values normalized to unit 2-norm and equal for all users. 

It is immediate that ^{G, r) is included in the feasible performance set (|32|) of the parallel 
channels and has the property that / G fi{G, r) implies af E ^{G, r), a < 1. 

The inequality (l58l) contains the tightest proposed bounds which utilize the extensions 
9'^, 9^ of the Lovasz function and Delsarte number. Since the intricacy of these bounds lies 
evidently in the structure of the sets ^{G, r), J^{G, r), we proceed by proving some loosened 
lower bounds which together with the looser lower bound (l30l) lead to our central insights. 

Corollary 3: Given any G = (/C, S) and r E we have 

max mm 7^ > max -. — \ , 

VV'=\,„a_r'^iB){wiwi ~B)+I 

with w such that f l27]) . with N E N satisfying X42\l . and, given A defined as in Corollary [7] 

HG,r) = {fEM!l:{fJ)<{fJ)^f = UPk),kElC, for some P eV{A)}. 

Moreover, given a particular ff3^ . this further implies 

. (akJkiPk))^ ^ /^(^^)max^g^(G,^)(/,/) 
max mm ^ > max — — . 

{A,P)eA{G,r)xV{A) kelC Yk VeR^^^: 9^{G,w) 

VV'=X„,ay-'^iB){w^w^^-B)+I 

Proof: The definition of TiG, r) implies that / E J-'{G, r) if and only if 

(/,/)<max5 subject to B{6) C HkeKif = fk{Pk) ■ P e V{A)}, (59) 

o>0 

where we define B{5) = {f E : (/, /) < 5}. On the other hand, we can write by the 
definition of ^{G, r) (quite redundantly) 

max (/,/)=min5 subject to B{5) ^ {f\k^K. {f = fk{Pk) ■ P ^ V{A)} 

n{/ = a/:«>0}), 

where / = argmin^gjgiv min^g^((; ^) jgjjiv max^e^ Now, as any / is arbitrarily 

nonnegatively scalable (that is, the latter set in the constraints in (l60l) is a ray in M.^), 
it is implied further by (l60l) that 

max (/, /)>min(5 subject to B{5) ^ B{5), /-/jix 

f(iP{G,r) &>0 > 



for any 5 satisfying the constraints in (|59l) . Thus, by (|59l ), (|6T| ) it follows finally 



(/,/)< max (/,/), /G^(G,r), 

/e^(G,r) 

which is, according to Proposition H sufficient for the proof for any given G = (/C, S) and 
r e M:^+. ■ 
By (|30l ) and Corollary [3] we have now 

'^^'^feHG,r){f^f) _ min^g^(g ,.)(/,/) 

o~/ \ — IXlclX ~ ~ 

e^{G,w) BeBHG,u,) \max[B) 

^ . (akJkiPk))^ . 
> max mm k > 

{A,P)eA{G,r)xV{A) keK 7^ (62) 

fi{RV)maXf^^^G,r){f,f) . ^RV) T^^^feHG,r){f , f) 

ITlcLX ; ; ^ 



VV'=X„,ax^^ {B)(wi wi ~B)+I 



with RV in the last expression as any nonnegative factor of the particular matrix (|57l) . with 
(|43] ). From the proof of the corollary it is evident that ^{G, r) can be interpreted as the 
largest, say, ball (in the 2-norm) of performance function values, equal for all users, included 
in each user dimension of the feasible QoS set (|32|) of the parallel channels. On the other 
hand, recall that the hull ^{G,r), determining the optimistic bound in (|62l ), contains any 
such user dimension of the feasible QoS set. Thus, the ball ^{G, r) and the hull ^{G, r) 
determine the interval (|62|) of candidate max-min fair performance values in terms of the 
structure of the feasible performance set of parallel channels; that is, in terms of the structure 
of the set of allowable power allocations 'P(A) and the features of the QoS functions f^, 
k E IC. In precise terms, the only such feature which is decisive for the bounds (|62l ) is the 
(squared 2-norm of the) minimum gap between T{G, r) and ^{G, r), measured as 

min {f,f)- max (/,/). 

Such gap is visualized, together with the hull !F{G, r) and the ball ^{G, r) for exemplary 
instance of parallel channels in Figs. |4] and [^. 

Consider now specifically the case of limitations of transmit powers at any time (in a 
frame), e.g. constrained transmit power of any user or constrained joint power budget of all 
users resulting in (l4al) . respectively. As in such case V{A) = V, A E A{r), it is readily 
seen that also ^{G, r) = ^{r) for any sharing graph G (i.e., T{G, r) is independent of 
the induced sharing graph). As a consequence, the influence of the features of the channel 
sharing policy on the interval (|62l ) of candidate values of max-min fair performance is 
in such case completely separated from the impact of the allowable power allocations. 
The combinatorial properties of the sharing graph G govern the inner bounds in (|62|) via 
the minimum achievable eigenvalues \,„ax{B) among matrices B E B'^{G,w) and the 
normalized eigenvalues X,„axiB) / fi(RV) among matrices B E B^{G, w); the normalization 
is by the (values of) the metrics /i of the associated factors of (|57] ). Analogously, the outer 
bound behavior is described by the function values 9'^{G, w) and normalized function values 
9^ {G , w) / fi(RV) , for the sharing graph G and the vector of squared user performance 
requirements per assigned channel fraction w, where the normalization is now by the metric 
1^1 of the corresponding factor of (1571 ) such that \nax{B) = 9^{G,w). Thus, the tightest 
pessimistic bound is obtained for a matrix B E B^{G,w) and a nonnegative factor RV 
of (l57l) which provide the minimum normalized eigenvalue \max{B) / ^{RV) . Similarly, the 

'Note here that the main results of this work are not proven to hold for the parallel channels instances from Figs. |4l |5] 
as the condition N > K is violated in these cases. Figs. [4] |5] serve, however, only as an exemplary visualization of the 
notions. 



fk2{pk2) = log(l + 




Fig. 4. The user dimensions of the feasible QoS set of parallel channels under the per-user power constraints (4b) 
and capacity (8) as performance function (left hand side), as well as the resulting hull ^{r), the ball ^{r) and the 
gap S = niinj,g^jj.j {/, /) — max^rg^p^^) (/, /) (right hand side). We simulated the parallel channels with K — 4 users 
fc = 1,2, 3, 4 accessing N — 2 channels n = 1,2. The channels hkn and the variances atn were picked randomly from 
uniform distributions resulting in an average Signal-to-Noise Ratio of 6 dB. 
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Fig. 5. The user dimensions of the feasible QoS set of parallel channels under the per-user power constraints (4b) and 
mean square reliability (7) as performance function (left hand side), as well as the resulting hull f{r), the ball T{r) and 
the gap 5 = iiiinj,g^jj.j {/, /) — maxjgjp(^) (/, /) (right hand side). We simulated the parallel channels as in Fig. 4, but 
for an average Signal-to-Noise Ratio of 9 dB. 



outer bounds in (|62l) are tightest for a factor RV of the particular (1571) . with (|43l) . which 
maximizes metric ^u. 

The outer bounds in (|62|) can be made in some sense symmetric whenever there exists a 
nonnegative factor RV of the matrix ( 1571 ) satisfying ( |43l ) which has row sums not exceeding 
r and each column sum no larger than 1 : In fact, as it is immediate from the Definition IH 
that then n{RV) > 1, we can embrace the max-min fair performance according to 

minj,^(^^,)(/,/) {akjkjp,))' max^g^(^,,)(/, /) 

-— -— r ^ max mm k ^ TrTTFi \ > \^J) 

where we have purely spectral dependence on the sharing graph G in the form of functions 
(O, (ITtI) . By the theory of matrix scaling [49], [48], the existence of such particular factor 
depends on the pattern of its zero entries, which is shown by the following paraphrased result 
from [53]. 

Proposition 5 ([53]): Let denote by RV{C\M.), with C C IC, M. C Af, the submatrix of 
RV G M;^^^ which is obtained by deleting all rows k E C and all columns n E AA from 
RV. Then, we have ji{RV) > 1 if 



^ffc<^c„ if RV{}C\C\M) = 0, RV{C\Ar\M) ^0, 
5^rfc=5^c„ if RV{IC\C\M) = 0, RV{C\Af\M) = 0, 



(64) 



holds for some (r, c) < (r, 1). Thus, the max-min fair performance satisfies ( |6?1 ) if there exists 
a nonnegative factor RV of the matrix ( 1371) . such that ( |4?1) and ^64\) for some (r, c) < (r, 1). 
Recall here from Lemma [21 that any matrix (1571) is a feasible matrix of some subgraph of the 
sharing graph. Thus, by the above proposition, the existence of a factor ensuring ii{RV) > 1, 
depends on the existence/nonexistence of certain edges in the sharing graph. 



C. Role of scalings 

The row and column sums of factors of the certain feasible matrix (1571) of some sharing 
subgraph influence the max-min fair performance in a specific way, which we show more 
explicitly here. Proposition |71 in Appendix O provides a technical alternative version of the 
bounds from Corollary [3l and we simplify it in the following. By the proof, one can readily 
see that the bounds from Proposition |71 are slightly tighter than those from Corollary [3l at the 
expense of higher complexit)0 Under apriori setting y = 1 in Proposition |71 and using the 
definition of w(x) and ( [T51 ), we obtain a more insightful, loosened version of the bounds: 
Together with (|30l) , we yield then precisely 

min^g^((5 ) (/, /) min_^g^(g ) (/, /) 

— max 



e^{G,w) BGBHC-w) Xmax{.B) 



\2 



. . {a-kJkiPk)? . 
> max mm ^ > 

{A,P)£A{G,r)xV{A) fce/C 7^ 

minfcg^Xfcmax^g^(c;,^)(/, /) minfegK;a;fcmax^g:^(G' )(/,/) 

UlclX - - - - 1 - 

B&BHG,w), Kax{B) ~ e^{G,w) 

VV'=\„a,-HB){'w^'w^' -B)+I 

(65) 

where in the last expression we can take any x satisfying {x,l) E X(RV,r,l) for any 
nonnegative factor RV of the matrix (1571) for (l43l) (note here that for ^{G, r, y), y E 



'"Note also that the proof of Proposition [7] allows for an even tighter bound formulation which generalizes Proposition 
|4j set T{G, r, y) has to be merely replaced by f(G, r, y) given in ( I82t . 



defined in Proposition |7] we have ^{G, r, 1) = T{G, r)). If the constraints of transmit 
power at any time in a frame are considered (e.g. either of (|4al)). then one can see the 
same separate impact of allowable power allocations and the channel sharing combinatorics 
on the bounds dM]) as in the case of ([621): In terms of V{A) = V, A e A{G,r), the 
interval of max-min fair performance values is determined by the distance between the 
largest included ball ^{G, r) and the hull T{G, r) of each user dimension of the feasible 
QoS set. Independently, the minimum achievable eigenvalues \nax{B) within sets (fT3] ) and 
(fTSi) . or the spectral characterizations 6*^ and 6^, govern the bounds in terms of the channel 
sharing topology expressed by the sharing graph G. By the definition of X{RV , r, 1), it is 
further evident that the outer lower bound in (|65] ) is a linear function of the minimum row 
scaling factor which is required to scale a nonnegative factor RV of (1571) down, until each 
column sum does not exceed unity and the vector of row sums is no larger than r. This 
leads to a conclusion that the outer bounds (l65l) embrace the max-min fair performance value 
as tightly as possible if such a nonnegative factor RV is taken which has componentwise 
smallest row sum vector relative to r. Note that, as (1571 ) represents some feasible matrix of 
some sharing subgraph G' C G, the row sums of a factor of (1571 ) are determined by the 
channel sharing combinatorics, that is, by the existence/nonexistence of certain edges in the 
sharing graph (and by the vector w of squared user performance requirements per assigned 
channel fraction). 

By the bounds (1651 ) it can be again seen that a symmetric embracing of max-min fair 
performance according to (|63l ) is implied whenever there exists a factor RV of (1571 ). for 
the particular (l43l) . which has all column sums no larger than unity and all row sums 
componentwise not exceeding r (see Proposition [51): In fact, in such case we can find a 
particular x such that min^gyc^^fc > 1 among all (cc, 1) G X{RV,r, 1). 

Consider now the complementary simplification of Proposition |71 where a; = 1 is set 
apriori. Then, together with (l30l ) we yield immediately 



max 



6'^{G,W) B&B-^{G,w) Xmax{B 



\2 



. {akJkiPk)y ^ 
> max mm 7^ > 

(A,P)£A(G,r)xV(A) fce/C 7^ 

T^^^f€HG,r,y)iy ° f^y ° f) > '^^^feHG,r,y){y ° f^y ° f) 



(66) 



where in the outer lower bound we can choose any y such that {y, 1) G X{RV, r, 1), with 
RV as any nonnegative factor of (1571) . where (l43l) . It is evident from the definition that the 
set ^{G, r, y) consists of performance function values which 

• are equal to, say, / G for all users accessing the parallel channels and, 

• when weighted by y in the sense y o f, are included in each user dimension of the 
feasible performance set of parallel channels. 

By analogy to ^{G, r), we can interpret the set ^{G, r, y) as a kind of largest ball which is 
included in each user dimension of set (l32l) . but which size is measured in a weighted (by y) 
Euclidean norm. Thus, the interval of max-min fair performance values (l66l) is influenced by 
the structure of the set V{A) and functions fk, k E JC, through the included weighted-norm 
ball r, y) and the hull T{G, r) of any user dimension of set (l32l) : the impact is purely 
via the weighted norm gap 



min (/,/)— max _ (/,/), 

/Gjf(G,r) {l,y)(iX{RV,r,l),f(G,T,y) 



where RV is a factor of (1571) for B E B^{G, w) achieving the tighter lower bound, or a 
factor of (1571) for (|43l) when the outer lower bound is considered. 

Recall that under constraints on transmit power at any time in a frame, such as (l4a1) . we 
have ^{G, r) = ^{r) and !F{G, r) = !F{r) regardless of G, and thus the power allocations 
and channel sharing graph influence separately the numerator and denominator of the bounds 
(l66l) . In addition to the impact of channel sharing combinatorics through the minimum of 
\nax{B) within (1T3l) and (1T51) (respectively, via the Lovasz function and Delsarte bound 
extensions 6*^, 6^), we see that the looser lower bound in (|66l) is proportional to the weighted 
squared 2-norm of y subject to (1, y) G X{RV , r, 1). Thus, the lower bound scales bilinearly 
with the vector of scaling factors which are needed in column-wise scaling of a nonnegative 
factor RV of (1571) to attain row sums and column sums componentwise not exceeding (r, 1). 
Again, as (1571) is a feasible matrix of a certain subgraph of the sharing graph, the column 
sums of RV are determined by the channel sharing topology and by the vector w. It can 
be observed that a factor RV which achieves smallest possible column sums is desired to 
provide as tight as possible outer interval of max-min fair performance values in (l66l) . The 
bounds (|66l) confirm the conclusion that we have the symmetric bounds (|63l ) whenever matrix 
(1571) satisfying (|43l) has a nonnegative factor with row sum vector no larger than r and no 
column sum exceeding unity (in this case (|63l ) is implied by (|66l ) by taking y = 1, and we 
also have f{G, r, 1) = f{G, r)). 

We can finally conclude that each of the inequalities proposed so far allows us to embrace 
by bounds also the nonrestricted max-min fair performance of parallel channels, i.e. the 
max-min fair performance when no sharing graph is given apriori. For instance, (|62l) implies 
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9'^{G,W) B&B-^{G,w) X,nax\B) 

> max mm ^ > 
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where G is a max-min fair sharing graph, i.e. a graph induced by the max-min fair sharing 
policy A, such that (that is, we have G = G{A) and A e A{G, r)). Clearly, the other 
bounds ( 1581 ). ( |63l ). (1651 ). ( |66l ) give rise to analogous enclosing of graph-nonrestricted max-min 
fair performance, when a max-min fair sharing graph is incorporated. 



D. Relation to the interference channel 

By (l66l) . one can recognize an interesting relation between max-min fair performance in 
parallel channels and the (weighted) throughput optimization in the interference channel 
considered in Section IV-AI Lemma HI and the definition of X{RV ,r,l) make evident 
that vector y in the lower bounds in (|66l ) corresponds to a certain power allocation in the 
associated interference channel. 

Corollary 4: Let an interference channel with user population M have an interference 
matrix {{RVy 0)' G M^^^, describing the interference among users according to [46], 
which corresponds to any nonnegative factor of d37l ) such that ^ |??1 ). Then, y in the outer 
lower bound is a power allocation in such interference channel which maximizes the 
weighted throughput 

k&K. 

with additive logarithmic power penalty terms {ck — fk)\ogZk, k E fC, and Ck\ogZk, k G 
J\f\]C, for some (r, c) < (r, 1). 



Precisely, by Lemma IH the weight vectors (r, c) in the throughput function collect row 
and column sums obtained under columnwise scaling of RV by y. The interesting point is 
that the throughput-optimal power allocation y in the described interference channel influ- 
ences the pessimistic bounds on max-min fair performance in the related parallel channels. 
For instance, the outer lower bound (l66l) becomes tighter if the taken nonnegative factor 
{(RVy 0)' of the matrix (1571) for the particular (|43l ) represents such an interference matrix 
of the associated interference channel, which enforces higher user powers for optimizing 
the weighted throughput from the corollary. Recall here that (l57l) is a feasible matrix of 
some sharing subgraph G' C G, so that the candidate interference matrices of the associated 
interference channel depend on the channel sharing topology in the original parallel channels. 

E. Role of sharing graph cycles 

What is apparent in all proposed inequalities enclosing the max-min fair performance 
so far, is the difference in the dependence on the channel sharing combinatorics between 
the upper and lower bounds. Upper bounds depend on the given sharing graph G (and 
weight vector w) through the minimum eigenvalue \nux{B) among matrices B E B'^{G, w), 
respectively through the associated value of the function 9"^. The lower bounds depend on 
the channel sharing policy via the minimum of A,„„,v(-B) among matrices B from the smaller 
set B'^{G, w), respectively via the value which the function 9^ assumes for G and w. By the 
recent results on completely positive graphs, we can, however, unify the dependence on the 
sharing graph for a large class of sharing graphs/topologies. 

Proposition 6: Let G = {IC,S) be any sharing graph with K < A or including no odd 
cycles longer than 4. Then, the bounds from Propositions |?] and Corollary \3\ and the 
bounds ^58\l . ^62\i . d6?l) . d65l) . d66l) are satisfied with 

B\G,w) = B^{G,w), and thus, 9^ {G,w) = 9^ (G, w) . 

The proposition is an immediate consequence of Lemma [2] and the definitions (fT4l) . (flTl) . 
The key to the above identity of 9'^ and 9^ is that, for any sharing graph G with no more than 
K = A vertices or no odd cycles longer than 4, any feasible matrix ( 1571) of a sharing subgraph, 
for any B E B'^{G,w), is completely positive and not only doubly nonnegative (see proof 
of Lemma 121). Proposition [6] implies that whenever the parallel channels are accessed by no 
more than K = A users, the value of the function 9'^ assumed for the sharing graph G (and 
vector w) is a sufficient characterization of the sharing policy for enclosing the max-min fair 
performance from above and from below, according to (1581 ). (|621) . (|63l) . (1651) or ( |66l ). Similarly, 
the value 9'^(G, w), for the given sharing graph G, is a sufficient description of the channel 
sharing for the proposed bounds (I58l) . (l62l) . (l63l) . (l65l) . (l66l) on max-min fair performance, 
when there is no odd chain of more than K = 4 users such that any two subsequent users 
share some channel and the last user shares a channel with the first one (this makes up a 
cycle in the sharing graph). In particular, we have such property when the users accessing 
the parallel channels can be partitioned into no more than M = 4 groups such that no pair 
of users within one group is allowed (or able) to share a channel; for instance, due to certain 
constraints on traffic class processing or hardware. The channel sharing is represented in 
such case by an A/-partite sharing graph, M = 2,3,4, with particular examples depicted in 
Fig. [3l Two parallel channel instances of this type, and thus such that the bounds (1581) . ( |62l ). 
(l63l) . (l65l) . (l66l) are determined solely be the function 9"^, were presented in Examples [6l and 
|71 In the multi-user multi-carrier channel from Example [6l certain user constellations are not 
allowed to share channels due to regulations on traffic processing, while in Example |71 the 
sharing of channels within some user classes is prevented/undesired because of excessive 
difference of delay times. 

To summarize, we note that whenever the parallel channels are shared according to 
any sharing graph G from Proposition [6l the proposed bounds enclosing the max-min fair 
performance are determined by the spectral properties of the channel sharing combinatorics 



via some value of X^axiB) among matrices B E B'^{G,w), respectively via the value of 
9'^ assumed by graph G and vector w collecting squared user performance requirements per 
assigned channel fraction. The structural features of the sharing topology have impact on 
the bounds through the metric ji, or row-sums, or column sums of a nonnegative factor of 
a feasible matrix (l57l) of some sharing subgraph. As far as transmit power constraints at 
any time (in a frame) are considered, e.g. ( |4al) . the impact of the (set of) allowable power 
allocations and the curvature of QoS functions is decoupled from the influence of the sharing 
graph; it is mirrored by the gap separating the hull ^{r) from the largest included ball ^(r) 
of each user dimension of the feasible QoS set. 

We close the discussion on the max-min fair performance by discussing the issue of 
the channel ensemble. It is evident from Corollary [T] and Propositions HI |7] that the max- 
min fair performance of parallel channels can be enclosed by bounds (|58] ). (|62l ), (|63] ). (|65] ), 
(l66l) whenever the number of accessed parallel channels satisfies (l42l) . This means that the 
proposed bounds apply to non-overloaded parallel channels for which the (cardinality of) 
channel population exceeds the (cardinality of) user population K at least by the factor 
{K + l)/2. From the proofs of Propositions HI |7] it is evident that such condition results 
from the use of the general nontight bound on cp-rank of matrix (157]) implied by Lemmas 
[H [2l As a consequence, the class of parallel channels instances satisfying (|42|) can be 
generalized, depending on the particular matrices B E B^{G, w) achieving the lower bounds 
in Propositions |4l U\ Precisely, the bounds (1581 ). (|621 ). (|63] ). (|651) . (I66l) apply, more generally, 
when 

N > max{K + 1, 4){\,„ax~^{B){w^w^' - B) + I)} 

is satisfied for the corresponding matrices B E B^{G,w) in the lower bounds. In other 
words, the proposed bounds apply, more generally, when the channel population exceeds the 
user population K by a factor no smaller than max{K + l,(f)}/K, with as the cp-rank of 
matrix (l57l) . for B achieving the lower bound of interest. 

VI. Characterization of some fair policies 

The proofs of the lower bounds from Propositions |4l|7] are constructive, that is, they contain 
implicit specifications of certain parallel channels policies. This allows us in this section 
to derive some algorithms for the computation of fair policies in the case of predefined 
topology, or equivalently graph, of parallel channels sharing. A fair policy is understood 
here as a policy which ensures user performance of any user be no worse than some 
specified pessimistic bound. According to Definition [H a predefined sharing graph means 
predetermined binary relations consisting in sharing/no sharing of channels by the single 
user pairs. We already explained in Section UlI] that the predetermination of channel sharing 
topology can be motivated by regulations on processing of different traffic classes, e.g. in the 
manner as in the multi-user multi-carrier channel from Example [6l The fixing of a channel 
sharing graph can be also necessary under certain constraints on hardware and/or signal 
processing, similarly to the Example |7] of parallel channels. 

In order to simplify the presentation, we assume that the predefined sharing topology 
results in a sharing graph with no odd cycles longer than 4, so that we have the equivalence 
from Proposition [6] throughout this section. Also, we restrict our attention to constraints on 
transmit power, e.g. by assuming individually constrained user power or constrained joint 
power budget of users at any time (in a frame) according to (l4al) : As a consequence, in what 
follows we have ^{G, r) = !F{r) and r) = ^{r) regardless of the sharing graph G. 
It is, however, easily verified that all the algorithmic concepts proposed in the following are 
straightforwardly extendable to the case of energy constraints (per frame). 



A. Fair policy as orthonormal-like representation 

Using the conventional optimization formulation, the problem of ensuring max-min fairness 
under given channel sharing topology can be written as 

{dkJkiPk)) / {A,P)eA{r)xV _ 

mm max , subject to < , ' , . , (67) 

{A,P) keK Ik \ {ak,ai)<Q, {kJ)^S, 

where the set £ is such that K?\8 collects all user pairs which are not allowed to share a 
channel or, equivalently (Definition [T]), G = {}C,S) is a given sharing graph0. Conventional 
optimization methods (e.g. interior point methods [54]) allow for a global solution basically in 
the case of convexity of the problem. Such feature is prevented in (l67l) since a bilinear form, 
used in the constraints, is not a convex function. Additionally, we consider very general 
performance functions fk, k E JC, and arbitrary constraints on transmit power, so that a 
standard method solution of (|67l) is expected, in general, to be only local. In this light, 
resorting to efficient computation methods of (suboptimal) fair policies seems to be an 
attractive alternative. 

One possible method is implied in the proof of Proposition |7] by the inequality (a refor- 
mulation of the first inequality in (fSTT )) 

{y° f,y° f) ^ . (afc,/)^ 

< mm 



[A, c) e A{r) X 



:) = argminmax — subject to I {a^^ai) < 0, {k,l) ^ S 

{y°f,y°f) [ {yoc,yoc)<l, 

(68) 

given any B G B'^{G,w{x)), with map z 1-^ w{z), z G defined in Proposition0 El 

Hereby, any vectors x, y satisfying (a:;, y) G X(RV, r, 1) for (l45l) can be chosen. It is evident 
by (fTT)) that the problem in (l68l) is closely related to the computation of an orthonormal 
representation and a unit vector which achieve the value of the Lovasz function (fTT)) (recall 
the definition of orthonormal representation from Section Ull]): In (|68l) . the unit vector c is, 
however, considered in weighted norm and is additionally restricted to be nonnegative, while 
the constraints on A are expressed in 1-norm. The complexity of the problem in (l68l) is 
significantly reduced in relation to the original problem (l67l) . We can restate this problem as 
an instance of so-called bilinear program by replacing the objective by some variable s and 
by adding the inequalities {ak, c)/7fc — s < 0, k E JC, to the constraints. Although a bilinear 
program does not represent a convex problem, there exists a variety of efficient methods for 
its global and local solution; without giving further details we refer for a selection of such 
methods to [55], [56] and references therein. 

Clearly, in the orthonormal-like representation (A, , =) obtained from the bilinear 

V {y°hy°f) 

program in (1681) vector / is arbitrarily scalable by a > 0. Due to our assumption a 
power allocation P E V satisfying 

fkiPk) = «/, kE /C, (69) 

always exists and is trivially constructed whenever « > is chosen sufficiently small: Under 
an appropriate a, any user k E JC accessing the parallel channels simply assigns on any 
channel n G A/" a transmit power pkn which achieves performance a/„ and the resulting 
power allocation P remains allowable. By iterative increasing of a in suitably small steps, 

"Note, that the inequahty in the second constraint in l |67t is equivalent to equality as nonnegativity is implicit from 
A G A{r). 

'^From the objective of the problem it is readily seen that the last inequality constraint can be replaced by equality. 



we achieve, with some accuracy, the particular largest a. for which (|69l ) is yet fulfilled for 
some P G P. Under a simple structure of the set of allowable power allocations, e.g. (l4al) . 
such value of a is often computable directly/non-iteratively once the performance functions 
fk, k E }C, are known. For such particular a we achieve the tightest lower bound in (|68l) 
among all a f inside the set JF(r, y) which is further smaller than the corresponding bound 
for any af G f{r, y) (recall dfcl), 

By (the proof of) Lemma HI a candidate vector y in (l68l) is computable as a solution 
of an unconstrained convex problem. As a first approach we prefer, however, to apply the 
simplification y = 1, which implicitly enforces x to satisfy (a;, 1) G X(RV,r,l). This 
results in the following simple procedure, for which r, 7^,, k E JC, and the set E of user pairs 
not allowed to share a channel (equivalently, sharing graph G = {IC,S)) are given as input 
parameters along with some suitably small a, 5 > 0. 

Algorithm 1: 

1: Find a sharing matrix A and vector f from ^6E\) . y = 1, by any bilinear programming 

method [55], [56]. 
2: Compute a power allocation P from d^Pl ). 
3: If P eV then set a ^ a + 6 and go to step 2, otherwise stop. 

With Proposition [6l the user performance of the obtained fair policy {A, P) is immediately 
evident from the proof of Proposition |7] (see bounds (1651)). 

Corollary 5: Given 8 such that G = {IC,S) has no odd cycles longer than 4, the policy 
{A, P) from Algorithm [7] satisfies the bounds from Proposition ^for y = 1, which implies 

e^iG,w) -fce^ 7| - e^{G,w) 

where {x, 1) G X{RV, r, 1) subject to PTl) and 

B = arg min A„,^(S). (70) 

Thus, min^gx: ^"'"■^''^^''-'^ is at most 

min_^g^(^)(/,/) -minfcg^Xfcmax_^g^(^)(/,/) 

e^{G,w) ' ^ ^ 

away from the max-min fair performance under given £. 

Fig. [6] provides an exemplary comparison between the user performance achieved by the 
policy from Algorithm [T] and the max-min fair performance. For the evaluated ensemble 
of parallel channels (with their sharing graphs) we observe a loss of about 20 % to the 
max-min fair performance. One can also show by simulation that such loss decreases if the 
differences between the user channel vectors hk, k E K,, and the differences between the 
variance ensembles cr^n, n E A/", of users k E K, diminish. In fact, this behavior can be 
recognized already from the feature (|69| ) of the policy from Algorithm [U which means that 
the resulting performance function value is the same for all users. Using Corollary [5l the 
same limit behavior can be also deduced from Figs. |4] and [5] since in the case of similar 
channel vectors and variance ensembles the (forms of) user dimensions of the feasible QoS 
set become similar as well and make the 2-norm gap between the hull ^{r) and the ball 
jF(r') vanish. As can be expected conversely, under variations between the channel vectors 
and variance ensembles of users becoming more severe, the loss of the policy from Algorithm 
[H increases. 

The advantageous complexity-performance trade off of Algorithm [1] becomes evident when 
the bilinear program in (l68l) and the original problem solution (l67l) are both computed by 
the same local optimization method. As shown in Fig. |7] for some selected parallel channel 
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Fig. 6. The comparison of user performance under the policy {A, P) from Algorithm 1 (dashed line) with the max-min fair 
performance (solid line), with sum-power constraint (4a) and with the capacity (8) as performance function. We simulated 
parallel channels with K — 4 users accessing N = 6 channels, equal user performance requirements 7^ — 1, k £ IC, and 
r = N/Kl. The sharing graphs G{j), 1 < j < 100, were picked randomly from all graphs with vertex set IC and edges 
occurring independently with probability 0.5. The channels hkn{j) and the variances Ukn{j), 1 < J < 100 were picked 
randomly from uniform distributions resulting in an average Signal-to-Noise Ratio of 20 dB. 



instances (and sharing graphs), the efficient and widely used Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) method may be attracted by highly suboptimal local optima of the original 
nonlinear problem. On the other hand, the values of the local optima of the bilinear program 
in (1681 ) are apparently much less scattered, so that the same BFGS method is able to find 
a good (local) solution (l68l) quite reliably. As a result, the worst user performance under 
policy from Algorithm [T] happens to be superior to the locally computed max-min fair policy 
(under given sharing graph). 

As a second approach to the scaling of RV, instead of the simplification y = Iwe can find 
a scalar scaling so that (1, yl) E X(RV, r, 1) for some y > 0. In this case, w{x) reduces to 
w defined as (|2TI) and thus matrix B, which gives rise to the nonnegative factorization (|45]) . 
needs to satisfy B E B'^{G, w). In the best case, a particular matrix (TTOI) is desired. Since the 
constraints determining the set (fT3l) are linear, the problem in (TTOl) corresponds to eigenvalue 
minimization over a polyhedron, which is a canonical problem in optimization theory and 
a variety of efficient solution methods exists [54]. For the nonnegative factorization of any 
given B E B'^{G, w), or the particular (iTOl) . we use one of the two celebrated methods which 
are proposed in [57] and are further extended and analyzed e.g. in [58], [59]. Precisely, for 
any m G N, we apply the particular form 

(V(m + 1)) - (^M).n V- (-^-M).JA„.r^(B)(^^^r-B) + J),, (kn)E]CxAf 

(72) 

of the factorization iteration from Theorem 2 in [57]. The sequence R'V{m), m E N, 
obtained by ([721) converges monotonically to a matrix RV which achieves a stationary point 
of the generalized KuUback-Leibler (KL) divergence between RVV'R and (|57]) (for the 
definition of this divergence and further discussion we refer to [57]). The minimization of 
the generalized KL divergence between a matrix and its factorization is an intricate problem 
with multiple local minima, so that RVV'R obtained from ([72] ) can happen to remain at 
a nonzero, but relatively small, generalized KL divergence to (|57]) . For this reason, we can 
resort also to alternative factorization iterations, such as the gradient descent method, which 
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Fig. 7. The comparison of user performance under tiie policy (A, P) from Algoritlim 1 using tlie BFGS method in step 
1 (dashed hne), and under pohcy {A, P) as a local BFGS solution to problem (67) (solid line), with sum-power constraint 
(4a) and with the capacity (8) as performance function. We simulated the parallel channels as in Fig. 6, but for K = 3, 
A*' = 4 and an average Signal-to-Noise Ratio of 10 dB. 



seem, however, to be inferior to the methods from [57] in terms of complexity-convergence 
trade off [57], [59]. 

The above discussion leads to the following second procedure which uses r, 7^, k G JC, 
and S as input data and some sufficiently small parameters a,6 > 0. 
Algorithm 2: 

1: Compute a matrix ( iTOl) by any convex eigenvalue minimization method [54]. 
2: Compute V by the iteration i[72\) . 

3: Compute the largest solution y > of the inequalities yl'RV < 1', yRVl < r. 

4: Find a sharing matrix A and vector f from ( |^g| ). y = yl, by any bilinear programming 

method [55], [56]. 
5: Compute a power allocation P from ( l^9l) . 
6: If P eV then set a ^ a + 5 and go to step 5, otherwise stop. 

By Proposition [6l Theorem 2 in [57] and the proof of Proposition |7] we have the following 
result on the user performance under the fair policy (A, P) from Algorithm [2l 

Corollary 6: Assume the generalized KL divergence between ( 1371) and RV, with V com- 
puted in step 3 of Algorithm^ be zero and let £ such that G = {fC,£) has no odd cycles 
longer than 4 be given. Then, the policy ( A, P) from Algorithm^ satisfies the bounds from 
Proposition ^for y = yl, with y computed in step 3, which implies 



Thus, min^gx: ^"'"'^''a^'''''^ is at most 



is at most 



(74) 



e^{G, w) 

away from the max-min fair performance under given £. 

Obviously, by adding more technicality. Corollary [6] can be extended to the case when the 
factorization iteration in step 3 happens to converge only locally, i.e., when the KL divergence 
between RVV'R and (|57l) does not vanish. 



The perfomiance and complexity-performance trade off of the policy from Algorithm |2] 
behaves, essentially, quite identically to the policy from Algorithm [2] (see Figs. [6l |7]). The 
potential nonzero KL divergence remaining after iteration by (|72l ) is hereby hardly visible. 



B. Fair policy from factorization 

To summarize so far, by the Algorithms [U [21 the solution of the original intricate problem 
(|67l) is replaced by some algebraic operations and the solution of canonical, more efficiently 
solvable optimization problems: The sharing matrix is obtained directly from the solution 
of a bilinear program, while the power allocation results from simple scaling (Algorithm 
[B, respectively, from the solution of eigenvalue minimization, nonnegative factorization and 
scaling (Algorithm O. As the price payed for this simplification, the resulting fair parallel 
channels policies are suboptimal, but achieve the worst user performance within the distances 
(ItTI) and (l74l) . respectively, from the optimum under given sharing graph G. 

The proof of Proposition |7] implies, however, that a bilinear program can be further 
exchanged here by nonnegative factorization and a solution of a simple equation system. 
The key step of the proof which gives rise to such alternative algorithm is the equality (see 

— r — -^7 — = 2 — ' keJC, B e B [G, w{z)), (75) 

for the given sharing graph G, for any B E B'^{G,w(z)), for map x t-^ w{x), x G 1^++, 
defined in Proposition |7] and for some z E A E A{G,r), f E related as follows. 

• Any sharing matrix A results from scaling of a factor RV, such that (l45l) . by a scaling 
(Z, Y) with {z, y) E X{RV, r, 1), where Zl = z,Yl = y. 

• Factor RV can be split as (|49|) . w = w(z), with X = (xi, . . . , xk)' as a factor of 
I - \max~^{B)B, where C = {c,. . . ^c)' E M^""^ is orthogonal according to (gV]) and 

determines / as , = = c. 

V {yof,y°f) 

Again, in the best case, a nonnegative factor RV of a particular matrix (TTOl ), w = w{z), 
obtained from canonical eigenvalue minimization [54], is desired. As above, such factor is 
computable by the version (|72|) of a factorization method from [57]. Once a factor RV is 
computed, vector / follows as a solution of a simple vector equation. Precisely, combining 

(|49l) with the definition , = = c and the orthogonality condition (|47]) shows that / is 

V {y°f,y°f) 
a solution to the equation 

y =-X„u,.^\B)wkz) = 0, fEM't, (76) 

V{y°f,y°f) 

where the definition VF2(z) = diag{W^ (z)), (W^ {z))^,^, = w^{z)f,, k E IC, is obvious. 

We are free to solve the equation (TTBl) by any available numerical method; we refer here 
to [60] for a wide selection of such methods. Any solution to ([761) is arbitrarily scalable 
by a positive a and from the discussion in Section IVI-AI it is clear how a power allocation 
P E V satisfying (l69l) is constructed for a sufficiently small a. Again, by gradual increasing 
the particular largest a is found, for which (l69l) yet holds for some allowable power allocation 
P E V. For such an a, a f achieves the value of the left hand side of (|75l) , which is further 
no smaller than the corresponding maximum value among all a f E ^{r, y). 

As a simplified approach to the scaling of RV , we find a scalar scaling which yields 
(1,?/1) E X[RV 1), for some y > 0. Since then w{z) reduces to w given by (12T]) and 
matrices B E B'^{G, w) have to be considered, the above discussion results in the following 
procedure (as above, S, r, and jk, k E IC, together with suitably small a, 5 > are given 
as input parameters). 

Algorithm 3: 
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Fig. 8. The comparison of user performance under thie policy {A, P) from Algorithm 3 (dashed line) with the max-min fair 
performance (solid line), with sum-power constraint (4a) and with the capacity (8) as performance function. We simulated 
the parallel channels as in Fig. 6, but for AT = 6 and N = 7. 



1: Compute a matrix ( |70l ) by any convex eigenvalue minimization method [54]. 
2: Compute V by the iteration I^TTt) . 

3: Compute the sharing matrix A = yRV, for the largest solution y > of the inequalities 

yl'RV < 1, yRVl< r. 
4: Compute a solution f to equation l^76\) by any numerical method [60]. 
5: Compute a power allocation P from ( 1^91 ). 
6: If P eV then set a ^ a + 5 and go to step 5, otherwise stop. 

According to Proposition [6l Theorem 2 in [57] and Proposition |71 the fair policy (A, P) 
computed by Algorithm [3] achieves the following user performance. 

Corollary 7: Assume the generalized KL divergence between ( 1371 ) and RV, with V com- 
puted in step 2 of Algorithm \3\ be zero and let £ such that G = {fC,£) has no odd cycles 
longer than 4 be given. Then, the policy (A, P) from Algorithm |2] satisfies the bounds 
from Proposition ^for y = yl, with y computed in step 3, which implies ( [73]) . Thus, 
min/cgx; ^"''"■^'4^''^^ is at most ^74^ away from the max-min fair performance under given 8. 
Fig. [8] shows an exemplary comparison of user performance achieved under the policy from 
Algorithm [3] and the max-min fair performance. It is evident that the average loss to the 
max-min fair performance is about 23 % for the simulated instances of parallel channels and 
their sharing graphs (thus, the potential nonzero KL divergence remaining after iteration (l72l) 
does hardly manifest itself in a gap to the performance of Algorithm [T]). By the feature (|69| ). 
or by Corollary |7] and the Figs. IH [5l we recognize again that such loss evolves analogously 
as in the case of Algorithms [T] and [2l it decreases with the user channel vectors and user 
noise variance ensembles converging to common values, and increases with the corresponding 
variations becoming stronger. 

VII. Conclusions 

This work allows for several novel conclusions on the behavior of the max-min fair perfor- 
mance in parallel channels, understood as the maximum attainable worst user performance. 
We assumed a very general performance function which is subject to the max-min fairness 
criterion; it includes the most celebrated functions in communications and information theory 
(capacity, spectral efficiency, decoder reliability) as very special cases. We succeeded in 



embracing the max-min fair performance by optimistic and pessimistic bounds which show, 
under constraints on transmit power, the same behavior as functions of the channel sharing 
topology. This lead to the first central conclusion that the max-min fair performance in parallel 
channels behaves as a special extension of the Lovasz function, or Delsarte bound, of a certain 
graph G (the sharing graph) characterizing the combinatorial topology of channel sharing 
among the users. An essential role is played hereby by the minimum spectral characterization 
\max{B) achievable within certain G-dependent sets B'^{G, w), B'^{G, w) with vector w as a 
parameterizing vector determined by the user performance requirements. When such spectral 
characterization of the channel sharing topology is obtained, the characterization of the real- 
valued subproblem of power allocation to users and shared channels by a simple 2-norm 
distance is sufficient for embracing the max-min fair performance by the proposed bounds: 
The influence of all properties of the allowable power allocations and all analytic features of 
the used QoS functions on the max-min fair performance is accumulated in a simple 2-norm 
gap between a certain hull and a certain included ball of the feasible QoS set of the parallel 
channels. 

Our results showed also that a key role is played by the existence/nonexistence of cycles in 
the sharing graph, interpretable as closed chains of users such that any two subsequent users 
in such a chain share a channel: We showed that under nonexistence of long odd chains 
of such type, the max-min fair performance is characterized by the minimum achievable 
\,nax{B) in the specific set B'^{G, w) as a function of the channel sharing topology (and the 
gap between the proposed bounds is equal precisely to the 2-norm gap between some hull 
and some included ball of the feasible QoS set). As a byproduct of our calculations, we 
also illustrated a relation of the max-min fair performance in parallel channel to the graph 
capacity and independence number of the graph describing the channel sharing topology. 

The constructive proofs of our bounds allowed further for the formulation of three novel 
power and time allocation algorithms for parallel channels with predefined channel sharing 
topologies (which is the case, e.g., under certain regulations/constraints on QoS class pro- 
cessing). The algorithms offer a nice performance-complexity trade off and incorporate some 
surprising techniques, such as nonnegative factorization. 



Appendix 

A. Doubly nonnegative and completely positive matrices 

Definition 5 ([44]): A matrix X G R^^^ is said to be doubly nonnegative, and we write 
X e ifDe M^^^ and D h 0. 

In simple words, a matrix X E D^^ is nonnegative in the conventional order > on M^^^ 
and in the partial order y on the set of symmetric matrices in M^^^. 

Definition 6 ([44]): A matrix X G R^^^^^ is said to be completely positive, and we write 
X G P^, if there exists some G N such that 

X = YY' for some Y G R;^^^. (77) 

The smallest number N for which we have is referred to as the cp-rank of X and is 
denoted as N = 0(X). 

Condition (fTTl) is frequently used in its equivalent form as 

N 

X = Y^ for some y,eR^, 1 <i<N, 

i=l 

where Y = (yi, . . . , yj^) is assumed. 

By Definition [6l it is readily seen that X G implies X G (P^' C D^^'). By the 
celebrated result from [44] it is further known that = D^^ whenever K < 4, while 
otherwise examples of matrices X G such that X ^ P^ can be constructed. 



B. Association schemes 



From the view of graph theory, the most accessible definition of an association scheme is 
based on the notion of edge coloring of a graph, as a partition of its edge set into vertex- 
disjoint edge classes. Precisely, an edge Af -coloring of a graph G = {}C,£) corresponds to 
the tuple (/C, {Sijfii), where {k, I), (m, n) G Si implies that k ^ m and / 7^ n [30]. 

Definition 7 ([61]): An association scheme with M associate classes on a set fC is an 
edge M-coloring of a (complete) graph G = (/C, K?) such that 

i.) for any 1 < k,l,m < M there exists so called intersection number p^] G N such that 

Pki = \{n e IC : (^i,n) e Sk, {n,j) e Si}\ whenever G £m., 

a.) for any 1 < k < M there exists qk E N such that qk = E Sk : i = n}\ for any 

n G /C. 

Hi.) Sk ^ 0, I < k < M. 

Definition 8: Given an association scheme (/C, {Si}fii) and any M. C {1, . . . , M}, we 
refer to C (Z IC as an M.-clique of the association scheme if i,j G C, i 7^ j, implies 

For any association scheme (/C, {Si}fii) and any its A^-clique C C IC, the unweighted 
Delsarte number can be formulated as the map 

{{IC,{S,}fi,),C,M)^ max l + (78) 

a£AH(IC,{S,}fi^),C,M)- 
l+(a,cr)>0 

where {{IC, {Si} fii), C, Ai) denotes the set of so-called inner distributions of the Ai- 
clique C and cr = (ai, . . . , (Ta/) G M^^ collects especially normalized eigenvalues of the 
adjacency matrices of graphs {}C,Sm), 1 <m < M, having a common eigenvector [39]. 



C. Additional bound formulations 

Corollary 8: Given N > K, any G = {IC, S) and any {A, P) G A{G, r) x V{A), 
r G we have 

{akJk{Pk))^ min^gKiv(/,/) 
min ^ ^^■'Y'''' < ^+ ^ = 0,1,2, (79) 

fce/c 72 e^{G{A),w{f)y 

where f w{f), f G M^, is such that, for any k E )C, 

Proof: The proof is a slight modification of the proof of Proposition 13] For any A G 

A{G,r), / G M^' and P G V{A) such that {ak, fk{Pk)) 7^ 0, k e JC (which by our 
assumptions on V{A) and fk, k E }C, exists), let Z = {zi, . . . , zk)' G M^^^, be given as 



^ _ Wk{f) n lk^{fJ) ^ , 

V (/, /) rk{ak, fk{Pk)) 

Then, 

lk\/wi{f){ak,f) li\/wk{f){ai,f) ^ lk^i{f,f){a'k,0'i) 
rk{ak, fk{Pk)) ri{ai, fi{pi)) rkri{ak, fk{Pk)){aiJi{Pi))' 

for any k,l E IC, so that by the definition of w in the case z = 0, 2 we yield again (l24l) . with 
w = w{f) and with {zk,Zk) = —Wk{f) + 'l!^^{!il'jllpj))'2 particular, while by Definition 
\T\ again (1251) for w = w{f) is satisfied. In the case of i = 1, the definition of w implies 
{zk, Zk) < -Wk{f) + ^|afe2(^)^ together with Definition H] also 



Zk, Zi) = ^J Wk{f)wi{f)- 



{zk, zi) < -VWk{f)wi{f), {k, I) ^S, k^l. 



Thus, given i = 0,2, we can write (|26l) with A = A and w = 'w(f) by the definition (fT3] ). 
while ^ ZZ' - (/, /)G(A) is satisfied for some B e B\G, w{f)) in the case i = 1. 
In either case (1271 ) is implied and up from dTT] ) the proof goes as the proof of Proposition 

m ■ 

Proposition 7: Given any G = (/C, S) and r G we /zave 

. {akJk{Pk)? ^ ^a^fenG,r,y){y°f^y°f) 

max mm ^ > max — , 

(a;,y)GA'(fiV,r,l), VgR^5^^: 



1 1 ' 



where x i— >■ w(x), x G M^i , is such that 



VV' = \„,ax~^{B)(w^W^ -B) + I 



where G N satisfies f l42]) . an J where, given A defined as in Corollary [7] we defined 
fiC, r,y) = {feR't:{yof,yof)<{yof,yof)^ 

f = fk{Pk),kelC, for some PeV{A)}, 
y G Moreover, given a particular 

B = arg min A„,,„(^), 

this further implies 

. {akJk{Pk)? ^ "^ax^g^(G )(yo/,7/o/) 
max mm 5 > max ^„ , „ — — . 

(Ai')e^(G,r)xP(yl) fcG/C 7^ (;r;,y)GA'(fiV,r,l), VeK^'x^: 9^[B,W[x)) 

VV' = \,„ax'^{B)(w^W^ -B) + I 

Proof: Under the substitution w = 'w{x), with an arbitrary x G the proof goes 

exactly as the proof of Proposition |4] up to the implication (|52|) for any / G such that 

— = c is satisfied for the particular c from (l50l) . Again, by Lemma [H it follows that we 

V (/./> 

can always find some {z,y) G A:'(i?V,r', 1) such that (|45|) . it; = 117(2;), and thus we have 
U G ^(G, r) for U = ZRVY with Zl = z, Yl = y. Writing now / = y o /, for some 
/ G M+, and U = {ui, . . . , uk)' and setting a; = 2, we have by the definition of the map 
X I— > w(x), X G that 



7I 7fc 7fc 



With (|52l). It; = it»(z), this yields 



ke)C. 



keJC, (80) 



for any / G M^^ such that , _ = = c is satisfied for c from (l50l) . By our assumption Q 

V {y°f,y°f) 

and the assumptions on Z^, A; G /C, we can always find a particular / such that / = fk{Pk)^ 
k e IC, for an arbitrary A G A{G, r) and some P G P(A), so that it is implied then with 
dSOl) that 

> mm max ^ — 

{A,f)eA{G,r)xRl k&K. {ak, j-} 

il{y°f,y°f) . il{y°f,y°f) 

= mm max ; — = mm max 

AeA{G,r),fm'l: k&IC {ak,fr AeA{G,r) kGtC (ttfc, /)2 

f=fk(Pk)MIC, for some PeV(A) 

(81) 



for any / G ^{G,r,y). Hereby, we defined (as a straight generalization of T{G,r)) 
:F{G, r, y) = {fe :f = fM, k e IC, for some P e V{A), 



. . {akJ? , (82) 

J = arg max _ max mm — — = 

A&A{G,v) keJC %{yof,yo f) 



Thus, it follows finally that 



\nax{,B) . 7I 

> mm max 



{y°f,y°f) {A,P)eA{G,r)xV{A) kdK. (a^, fkiPk))^' 

for any B e B^{G,w{z)), any V such that (l45l) . w = w{z), any y satisfying {z,y) e 
X{RV,r, 1) and > max{fs: + 1, (p{\„,ax~\B){w^z)w^z)' - B) + I)}. Hereby, by 
Lemmas [B [21 the last condition is implied by > K{K + l)/2 and additionally, along 
exactly the same lines as in the proof of Corollary [3l it is readily shown that 

{y°f,y°f)< max {yof,yof), fef{G,r,y). (83) 

This completes the proof of the first inequality of the proposition, while the second in- 
equality is obtained by taking a particular B with X,nax{B) = mai^B&B^iG w{z))^max{B) = 

e%G,w{z)). ' m 



References 

[1] L. Tassiulas and S. Sarkar, "Maxmin Fair Scheduling in Wireless Networks," Proc. IEEE Conference on Computer 

Communications (Infocom), Jun 2002. 
[2] B. Radunovic and J.-Y. L. Boudec, "Why Max-min Fairness Is Not Suitable For Multi-Hop Wireless Networks," 

Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Report, 2000. 
[3] D. K. Foley, "Resource Allocation and the Public Sector," Yale Economic Essays, vol. 7, 1967. 
[4] H. R. Varian, "Equity, Envy, and Efficiency," Journal of Economic Theory, vol. 9, pp. 63-91, 1974. 
[5] H. Boche, M. Wiczanowski, and S. Stanczak, "Unifying View on Min-Max Fairness, Max-Min Fairness, and Utility 

Optimization in Cellular Networks," Eurasip Journal on Wireless Communications and Networking, vol. 2007, p. ID 

34869, 2007. 

[6] A. Lozano, A. M. Tulino, and S. Verdu, "Optimum Power Allocation for Parallel Gaussian Channels With Arbitrary 

Input Distributions," IEEE Transactions on Information Theory, vol. 52(7), pp. 3033-3051, Jul 2006. 
[7] D. P. Palomar and J. R. Fonollosa, "Practical algorithms for a family of waterfilling solutions," IEEE Transactions 

on Signal Processing, vol. 53(2), pp. 686-695, Feb 2005. 
[8] G. D. Forney, Jr. and G. Ungerboeck, "Modulation and Coding for Linear Gaussian Channels," IEEE Transactions 

on Information Theor)', vol. 44(6), pp. 2384-2415, Oct 1998. 
[9] J. A. C. Bingham, P. S. Chow, and J. M. Cioffi, "A practical discrete multitone transceiver loading algorithm for data 

transmission over spectrally shaped channels," IEEE Transactions on Communications, vol. 43(2), pp. 773-775, Feb 

1995. 

[10] A. Lozano, A. M. Tulino, and S. Verdu, "Mercury/waterfilling for fixed wireless OFDM systems," in Proc. IEEE 

Radio and Wireless Symp. (RWS'06), San Diego, CA, Jan 2006, pp. 211-214. 
[11] I. Toufik and R. Knopp, "Multiuser Channel Allocation Algorithms Achieving Hard Fairness," in Proc. IEEE Global 

Telecommunications Conference (Globecom), 2004. 
[12] M. Wiczanowski, H. Boche, and S. Stanczak, "Characterization of Optimal Resource Assignments in the Framework 

of Blocking System Theory," in Proc. IEEE International Symposium on Information Theory and its Applications 

(ISITA), Seoul, Korea, Nov 2006. 
[13] , "Power Allocation and Resource Assignment in the View of Blocking and Antiblocking Polyhedra," in Proc. 

IEEE Information Theory Workshop (ITW), Chengdou, China, Oct 2006. 
[14] Y. J. Zhang and K. B. Lataief, "Adaptive Resource Allocation for Multiaccess MIMO/OFDM Systems with Matched 

Filtering," IEEE Transactions on Communications, vol. 53(11), Nov 2005. 
[15] D. R. Fulkerson, Blocking Polyhedra, ser. Graph theory and its applications. Academic Press, 1970, pp. 93-112. 

[16] , "Antiblocking Polyhedra," Journal on Combinatorial Theory, Series B, vol. 12, pp. 50-71, 1972. 

[17] L. Lovasz, "On the Shannon Capacity of a Graph," IEEE Transactions on Information Theory, vol. 25(1), pp. 1-7, 

1979. 

[18] M. Grotschel, L. Lovasz, and A. Shrijver, "Relaxations of Vertex Packing," Journal of Combinatorial Theory, Series 

B, vol. 40, pp. 330-343, 1986. 
[19] J. Korner and A. Orlitsky, "Zero-error information theory," IEEE Transactions on Information Theory, vol. 44(6), pp. 

2207-2229, Oct 1998. 

[20] X. Qin and R. Berry, "Distributed power allocation and scheduling for parallel channel wireless networks," in Proc. 
International Symposium on Modeling and Optimization in Mobile, Ad Hoc, and Wireless Networks (WIOPT), 2005, 
pp. 77-85. 



[21] C. Y. Wong, R. S. Cheng, K. B. Letaief, and R. D. Murch, "Multiuser OFDM with Adaptive Subcarrier, Bit, and 
Power Allocation," IEEE Journal on Selected Areas in Communications, vol. 17(10), pp. 1747-1758, Oct 1999. 

[22] H. Bolcskei, "MIMO-OFDM Wireless Systems: Basics, Perspectives, and Challenges," IEEE Wireless Communica- 
tions, pp. 31-37, Aug 2006. 

[23] 3rd Generation Partnership Project (3GPP), "Evolved Universal Terrestrial Radio Access (E-UTRA) and Evolved 
Universal Terrestrial Radio Access Network (E-UTRAN) (Release 8)," 3GPP Support Office, Valbonne, France, 
Technical Specification, 2007. 

[24] , "Evolved Universal Terrestrial Radio Access (E-UTRA); Base Station (BS) radio transmission and reception 

(release 8)," 3GPP Support Office, Valbonne, France, Technical Specification, 2008. 
[25] T. M. Cover, "Comments on Broadcast Channels," IEEE Transactions on Information Theory, vol. 44(6), pp. 2524- 

2530, Oct 1998. 

[26] N. Jindal and A. Goldsmith, "Capacity and Optimal Power Allocation for Fading Broadcast Channels With Minimum 

Rates," IEEE Transactions on Information Theory, vol. 49(11), pp. 2895-2909, Nov 2003. 
[27] D. N. Tse, "Optimal Power Allocation over Parallel Gaussian Broadcast Channels," in Proc. IEEE International 

Symposium on Information Theory (ISIT), 1997, p. 27. 
[28] D. Tse and P. Viswanath, Fundamentals of Wireless Communication. Cambridge: Cambridge University Press, 2005. 
[29] S. Verdu, Multiuser Detection. New York: Cambridge University Press, 1998. 
[30] B. Bollobas, Modern Graph Theory. New York: Springer- Verlag, 1998. 
[31] D. E. Knuth, "The Sandwich Theorem," Electronic Journal of Combinatorics, vol. 1, 1994. 
[32] D. M. Cvetkovic, M. Doob, and H. Sachs, Spectra of Graphs. New York: Academic Press, 1980. 
[33] J. H. Drew, C. R. Johnson, and R. Loewy, "Completely Positive Matrices Associated with m-Matrices," Linear and 

Multilinear Algebra, vol. 37, pp. 303-310, 1994. 
[34] N. Kogan and A. Berman, "Characterization of completely positive graphs," Discrete Mathematics, vol. 114, pp. 

297-304, 1993. 

[35] A. Berman and D. Hershkowitz, "Combinatorial Results on Completely Positive Matrices," Linear Algebra and its 

Applications, vol. 95, pp. 111-125, 1987. 
[36] C. Xu, "Nearly Completely Positive Graphs," Applicable Algebra in Engineering Communication and Computing, 

vol. 13, pp. 1-8, 2002. 

[37] C. J. Luz and A. Schrijver, "A Convex Quadratic Characterization of the Lovasz Theta Number," SIAM Journal on 

Discrete Mathematics, vol. 19(2), pp. 382-387, 1979. 
[38] A. Schrijver, "A Comparison of the Delsarte and Lovasz Bounds," IEEE Transactions on Information Theory, vol. 

25(4), pp. 425^29, Jul 1979. 

[39] P. Delsarte, "An algebraic approach to the association schemes in coding theory," Philips Res. Reps. SuppL, vol. 10, 
1973. 

[40] F. Chung and R. M. Richardson, Weighted Laplacians and the Sigma Function of a Graph, ser. Quantum Graphs And 

Their Applications - Contemporary Mathematics. American Mathematical Society, 2006, ch. 7, pp. 93-108. 
[41] A. Galtman, "Spectral Characterizations of the Lovasz Number and the Delsarte Number of a Graph," Journal of 

Algebraic Combinatorics, vol. 12, pp. 131-143, 2000. 
[42] R. J. McEliece, E. R. Rodemilch, and H. C. Rumsey Jr., "The Lovasz Bound and Some Generalizations," Journal of 

Combinatorics Information & System Sciences, vol. 3(3), pp. 134-152, 1978. 
[43] C. R. Johnson and M. Newman, "A note on Cospectral Graphs," Journal of Combinatorial Theory, Series B, vol. 28, 

pp. 96-103, 1980. 

[44] L. J. Gray and D. G. Wilson, "Nonnegative factorization of positive semidefinite nonnegative matrices," Linear Algebra 

and its Applications, vol. 31, pp. 119-127, 1980. 
[45] A. Berman, Completely positive graphs, ser Quantum Graphs And Their Applications - Contemporary Mathematics. 

American Mathematical Society, 2006, ch. 7, pp. 93-108. 
[46] S. Stanczak, M. Wiczanowski, and H. Boche, Resource Allocation in Wireless Networks. New York: Lecture Notes 

in Computer Science, Springer- Verlag, 2006. 
[47] A. W. Marshall and 1. Olkin, "Scaling of matrices to achieve specified row and column sums," Numerische Mathematik, 

vol. 12, pp. 83-90, 1968. 

[48] U. G. Rothblum and H. Schneider, "Scalings of Matrices Which Have Prescribed Row Sums and Column Sums via 

Optimization," Linear Algebra and its Applications, vol. 114, pp. 737-765, 1989. 
[49] D. London, "On matrices with a doubly stochastic pattern," Journal of Mathematical Analysis and Applications, 

vol. 34, pp. 648-652, 1971. 

[50] M. Schubert and H. Boche, QoS-Based Resource Allocation and Transceiver Optimization, ser. Foundations and 
Trends in Communications and Information Theory. Now Publishers, 2005, vol. 2(6). 

[51] U. G. Rothblum and H. Schneider, "Characterization of Optimal Scalings of Matrices," Mathematical Programming, 
vol. 19, pp. 121-136, 1980. 

[52] F. Murnaghan, The Unitary and Rotation Groups. Washington: Spartan Books, 1962. 

[53] R. A. Brualdi, "Convex sets of nonnegative matrices," Canadian Journal of Mathematics, vol. 20, pp. 144-157, 1968. 
[54] D. P. Bertsekas, Nonlinear Programming. Athens: Athena Scientific, 1999. 

[55] D. J. White, "A linear programming approach to solving bilinear programmes," Mathematical Programming, vol. 56, 
pp. 45-50, 1992. 

[56] H. D. Sherali and C. M. Shetty, "A finitely convergent algorithm for bilinear programming problems using polar cuts 

and disjunctive face cuts," Mathematical Programming, vol. 19, pp. 14-31, 1980. 
[57] D. D. Lee and H. S. Seung, "Algorithms for non-negative matrix factorization," Advances in Neural Information 

Processing, vol. 13, pp. 556-562, 2000. 



[58] P. O. Hoyer, "Non-negative Matrix Factorization with Sparseness Constraints," Journal of Machine Learning Research, 
vol. 5, pp. 1457-1469, 2004. 

[59] C.-J. Lin, "On the Convergence of Multiplicative Update Algorithms for Nonnegative Matrix Factorization," IEEE 

Transactions on Neural Networks, vol. 18(6), pp. 1589-1596, 2007. 
[60] J. M. Ortega and W. C. Rheinboldt, Iterative Solution of Nonlinear Equations in Several Variables. New York: 

Academic Press, 1970. 

[61] R. A. Bailey, Association Schemes: Designed experiments, Algebra and Combinatorics. Cambridge: Cambridge 
University Press, 2004. 



